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1 Introduction 



It is now forty years that Nicola Cabibbo formulated the notion of flavour mixing in the charged hadronic 
weak interactions^. The Cabibbo theory provides a consistent description of the muonic decay fi — > ev e v Ul 
the neutron /3-decay n — > pev e , and the strangeness changing transitions, such as the decays and the 
hyperon decays, in terms of a universal Fermi coupling constant Gf and a mixing angle, the Cabibbo 
angle 9c- Thanks to dedicated experiments carried out well over four decades, and impressive theoretical 
progress, in particular in the technology of the electroweak radiative corrections, we now have precise 
values for these fundamental parameters of nature •=* 

Gf = 1.16639(1) x 1(T 5 GeV~ 2 , 

(1) 

9c = 12.69(15)°. 

The Cabibbo theory ^ describes in the quark language charged weak transitions d — > u and s — > u 
involving the three lightest quarks u, d, and s. However, it was not able to account for the flavour 
changing neutral current (FCNC) transition s — > d, such as the K° - K° mass difference AMk- This 
outstanding problem with the Cabibbo theory was solved by the GIM mechanism ^ which required a 
fourth quark - the charm (c) quark. The GIM-construction banished the FCNC transitions from the 
tree level, relegating them to loops (induced quantum effects) where they found their natural abode. 
Thus, in the Cabibbo-GIM theory, AMk, as well as a number of |A5| = 1, AQ = transitions, such as 
Kl — > an d Kf — > 77, are well-accounted for in terms of Gf and 9c, and the mass of the charm 

quark m c l^l found later to be in the right ball-park through the discovery of the J/ip,tl>', ... resonances 
and charmed hadrons. Along with the GIM mechanism came also a (2 x 2) quark flavour mixing matrix 
characterized by the Cabibbo angle 9c- 

The final act in quark mixing came through the seminal work of Kobayashi and Maskawa (KM) ^ 
who enlarged the Cabibbo-GIM (2 x 2) quark mixing matrix to a (3 x 3) matrix by adding another doublet 
of heavier quarks (t, b). This matrix which relates the quarks in the weak interaction basis (d! , s', b') and 
the quark mass eigenstates (d, s, b), 

(d',s',b') = V C km {d,s,b) , (2) 
is called the Cabibbo-Kobayashi-Maskawa matrix VckMj an d symbolically written as 



Vud 


v us 


Vub 


Vcd 


Vcs 


Vcb 


Vtd 


Vts 


V tb 



VCKM = Vcd Vcs Vcb ■ (3) 

V Vtd Vts Vtb J 

Vckm is a unitary matrix, characterized by three independent rotation angles and a complex phase. 
The KM theory was formulated to incorporate the CP violation observed in Kaon decays in 1964 by 
Christenson et al. ®. In this theory CP symmetry is broken at the Lagrangian level in the charged 
current weak interactions and no where else. In principle, all the elements of the matrix Vckm are 
complex. In practice, only two of the matrix elements have measurable phases. But, this is sufficient 
to anticipate CP violation in a large number of processes, some of which are now being measured with 
ever-increasing precision in the K- and B-meson decays. 

In these lectures, I will summarize the current status of the magnitude of all the nine matrix el- 
ements \Vij\ of the CKM matrix and the weak phases entering in these matrix elements. To discuss 
this, a parametrization of the CKM matrix is needed. It has become customary to discuss the CKM 
phenomenology using the Wolfenstein parametrization^ 

/ 1-iA 2 A AX 3 (p-iri)\ 

Vckm * -A(l + iA 2 \ 4 v ) 1 - 5A 2 AX 2 , (4) 

\A\ 3 (1 - p - irj) -AX 2 (l + i\ 2 r]) 1 / 

where the four independent parameters are: A, X = singes P an d rj, of which r\ is what makes this matrix 
complex and leads to CP Violation. Anticipating precise data, a perturbatively improved Wolfenstein 
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parametrizationElwith p = p(l — A 2 /2) and f) = r](l — A 2 /2) will be used. This rescaling effects mainly the 
matrix elements Vtd, which now has the definition Vtd = AA 3 (1 — p— iff), and V u b — AA 3 (l + A 2 /2) (p — iff), 
and the other matrix elements remain essentially unchanged. 

As we shall see, a quantitative determination of these matrix elements requires, apart from dedi- 
cated experiments, reliable theoretical tools in the theory of strong interactions (QCD). These include, 
apart from the QCD-motivated quark models, chiral perturbation theory, QCD sum rules, Heavy Quark 
Effective Theory (HQET) and Lattice QCD, combined with perturbative QCD. To illustrate their im- 
pact, I will discuss some representative applications where a particular technique is the main theoretical 
workhorse. Further detai ls an d in-depth discussions can be found, for example, in the proceedings of the 
CERN-CKM workshops EE3 

These lecture notes are organized as follows: In section 2, I review the status of \V u d\ and |V US | and 
the resulting test of unitarity of the CKM matrix. In section 3, current status of \V c d\ and \ V ca \ is briefly 
reviewed. Section 4 describes in considerable detail the current measurements of the matrix elements 
\V c b\ and \V u b\ and the theoretical techniques used in arriving at the results. Status of the third row of 
Vckm is reviewed in Section 5. Section 6 summarizes the current knowledge of \Vij \ and the Wolfenstein 
parameters, including the phases of the unitarity triangle(s). This section also reports on the results of 
a global fit of the CKM parameters using the CKM unitarity and knowledge of the various experimental 
and theoretical quantities. CP violation in _B-meson decays is discussed in Section 7, and we restrict 
ourselves to the discussion of only the currently available results from the two B-factory experiments. We 
conclude with a summary of the main results and some remarks in Section 8. 



2 Current Status of \V u d\ and \V U8 \ 
2.1 Status of \V ud \ 

We start with the discussion of the matrix element \V u d\- The superallowed + — > + Fermi transitions 
(SFT) have been measured so far in nine nuclei ( 10 C, 14 0, 26 Al, 34 CI, 38 K, 42 Sc, 46 V, 50 Mn, 54 Co), sum- 
marized by Towner and Hardy ^5 As only the vector current contributes in the nuclear hadronic matrix 
element (p/; + \u'j f j i d\pi] + ), and the transitions involve members of a given isotriplct, the conserved 
vector current hypothesis helps greatly in reducing the hadronic uncertainties. Radiative (<5r and Ar) 
and isospin breaking (5c) corrections have been calculated. Of these 5c and 5r are nucleus-dependent. 
As the ft values (/ is the nuclear-dependent phase space and t is the lifetime) of the nuclear transitions 
are also nucleus-dependent, one usually absorbs the nucleus-dependent radiative corrections by defining 
another quantity 

Tt = ft(l + 5 n ){l-5 c ). (5) 

The Tt values mea sured in the nine transitions are indeed consistent with each other, with their average 
having a value 

ED 

Tt = 3072.3(9) sec [x7 d -c-f- = 1.1] . (6) 

The nucleus-independent radiative cor rection A# incorporates the short-distance contribution and has 
been calculated by Marciano and Sirlin 11^1121 The value of An depends on a parameter vtla which enters 
in the process of matching the short-distance and long-distance contributions. Taking vcia in the range 
ttiai/2 < niA < 2m ai, where tuai is the l" 1 meson mass, one estimates: Ar = (2.40 ± 0.08)%. With 
the precise determination of Tt and A^, the matrix element \V u d\ is derived from the expression: 

where K is a phase space factor, K — 2ir 3 ]xi2/m\, with m e being the electron mass, an d it s value is 
known to a high accuracy: K = (8120.271 ± 0.012) x 10~ 10 GeV- 4 sec. This yields a value El 

0.9740 ±0.0005. (8) 
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The error shown here is dominated by theoretical errors, contributed mainly by the (somewhat arbitrary) 
ch oice of the low energy cutoff in estimating A ^ and the nuclear-dependent isospin-breaking corrections 
S c HUH. The value listed by the PDG El from this method is: \V ud \ = 0.9740 ± 0.0005 ± 0.0005, where 
the added error reflects the PDG concern about the systematic uncertainty due to the nucleus-dependent 
radiative corrections. 

The other precise method of determining the matrix element \V ud \ is through the polarized neutron 
beta-decay (n — > pe~V^). The currently attained precision owes itself to the enormous progress made in 
having highly polarized cold neutron beams. For example, for the cold neutron beam at the High Flux 
Reactor at the Institut-Laue-Langevin, Grenoble, the de gre e of neutron polarization has been measured to 
be P = 98.9(3)% over the full cross-section of the beam Also, the neutron lifetime, r„ = (885.7 ±0.8) 
sec El is now measured to an accuracy of one part in a thousand. The charged weak current has a V — A 
structure and the hadronic matrix element can be parametrized as: (p | ii^p, ( 1 — 7s)c?|n) — Up'y^lgy + 
gAlb)u ni requiring the knowledge of gy and gA- (We have neglected a small weak magnetism contribution 
oc (pip — /i n )/(2m p )cr Mt/ g !y , where /i p and fi n are the anomalous magnetic moments of the proton and 
neutron, respectively, with q v and m p being the momentum transfer and the proton mass, respectively). 
However, in the neutron beta-decay, radiative corrections are under better theoretical control. 

The theoretical expression for determining \V ud \ from the neutron lifetime is: 

|Kd|2 = Cr„(l + 3*2)/* (1 + A R ) ' (9) 

where C = Gj^ml / (2ir 3 ) , x = gA/gv, / R = 1.71482(15) is the phase space factor including model- 
dependent radiative corrections'^, and the model-independent radiative correction A# has been specified 
earlier. The high accuracy on / R owes itself to the Ademollo-Gatto theorem, which makes the departur e 
of gy from the symmetry limit tiny, with current estimates yielding Sgy = 1 — gy — O(10~ 5 ) 

To extract \V ud \ from the neutron lifetime, one has to know gA/gv- This can be determined from 
the electron (/3) asymmetry or the e-D correlation in the decay of a polarised neutron. For example, the 
probability that an electron is emitted with an angle 9 with respect to the neutron spin polarization, 
denoted here by P — (a z ), is 

W{6) = 1 + (3PAcos9 , (10) 
where (3 is the electron velocity and the coefficient A depends on x = gA/gv'- 

A=-2^±4. (11) 
1 + 3x 2 y ' 

It is understood here that a small correction due to the weak magnetism has been included in extracting 
gA/gv from A. Thus, the measurements of r„ and A determine both gA/gv and \V u d\- However, 
currently the two most pre cise measurements of this quantity, namely gA/gv = — 1.2594 ± 0.0038^^1 and 
gA/gv = -1.2739 ±0.0019 ESI differ by more than 3<r, and hence the experimental spread in the values 
of gA/gv is currently the main un cert ainty in the determination of \V u d\ from the neutron /3-decay. Of 
these, the PERKEOII experiment ^1 yields a value \V ud \ = 0.9724 ± 0.0013, which, on using the PDG 
values for \V US \ (= 0.2196 ± 0.0026) and \V ub \ (= (3.6 ± 0.7) x 10~ 3 ) leads to 

\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 1 - Ax = 0.9924 ± 0.0028 . (12) 

The value Ai = (7.6 ± 2.8) x 10~ 3 differs from zero (the unitarity value) by 2.7<r. However, follow- 
ing t he advic e of the PDG and restricting to the experiments using neu tron polarization of more than 
90% | 15 | 21 | 22 | a recen t compilation of the experimental results yields ESI 

gA/gv = -1.2720 ± 0.0022, 

(13) 

|Kd|n decay = 0.9731 ± 0.0015 . 
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The result for \V u d\ from the neutron /3-decay is less precise than the one in ||SJ|, obtained from the 
+ — * + SFTs, though the two values of | | are completely consistent with each other. To improve 
the precision on \V u d\ from the neutron beta-decay, it is imperative to resolve the inconsistencies in the 
current measurements of gA/gv and determine this ratio more accurately. 

The third method for determining \V u d\ is through the 7r e 3 decay: tt + — ► ir°e + v e - This decay is 
governed by the vector pion form factor /7 w (t), where t is the transfer momentum squared. In the 
isospin limit. fj ±7r (0) = 1. An updated analysis of the radiative corrections to the pionic beta-decay 
has been recently undertaken in an elegant paper by Cirigliano et al. M including all electromagnet 
corrections of order e 2 p 2 (here e 2 p 2 implies both corrections of order e 2 p 2 and of order (m u — m,d)p 2 ), 
using the framework of chiral perturbation theory with virtual photons and lcptons. Accounti ng f or the 
isospin-breaking and radiative corrections, \V u d\ can be obtained from the following expression l«l 

|Kd|2 = WP%^M' (14) 

where 

K = G%M^ ± /(64w s ), 

I„(X t ,a) =/ w (A t ,0) {l + AI n (X u a)) , (15) 
/T ±7r °(0) = (l + ^ £/(2) )(l + 5 e % 3 ). 

Here X t is the slope parameter in the parametrization of the form factor f+ ±7T °{t) = f+ ±7l ° (0)(1 + 
X t t/M 2 + 0(X 2 )), I n (X t ,a) is a slope-dependent phase space integral, and S cw is what was earlier called 



Ar. Estimates of the various quantities in these expressions are^J 

°SU(2) 1U i 

^ = (0.46 ±0.05)%, (16) 
AI w {X t ,a) = 0.1%. 

The precision on \V u d\ is dominated by the precision on the quantity T^ e3 . , (i.e., the branching ratio 
BR(7r ± — * TT°e ± iy e [+j]). The present preliminary result of the PIBETA collaboration 1^51 is: 

BR(tt ± -> 7r°e ± i/ e ) = (1.044 ± 0.007(stat.) ± 0.009(syst.)) x 1(T 8 , (17) 

which is significantly more precise than the earlier most accurate measurement by Mc Farla ne et al. 
BR(?r ± -> 7r°e ± ^ e ) = (1.026 ± 0.030) x 10~ 8 . The PIBETA measurement yields a valued 

\VvdU* = 0.9771 ±0.0056. (18) 

This measurement of \V u d\ is almost an order of magnitude less precise than \V u d\ determined from the 
nuclear + — > + SFTs. One expects a factor three improvement in the value of \V u d\ at the end of the 
PIBETA experiment. 

Taken the three determinations of |V^<j| discussed here, the current world average of this quantity 

|Kd|wA = 0.9739 ± 0.0005 , (19) 

is essentially the same as in JHJ from the nuclear + — > + transitions. The impact of the neutron beta 
decay experiments on jV^I can be significantly enhanced if the experimental spread in gA/gv is resolved, 
and the resulting accuracy on this quantity improved. 
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2.2 Status of \V US \ 

The determinatio n of th e matrix element \V US \ from Kg 3 decays (with i = e,/j.) has been extensively 
reviewed recently l ^ l ^ l to which we refer for further details. The value for \V US \ quoted by the PDG 
is based essentially on the theoretical analysis of Leutwyler and Roos ^Hl done some twenty years ago, 
which yields: \V US \ = 0.2196 ± 0.0023. During the last couple of years, new analytical calculations of the 
radiative corrections have been reported by Cirigliano et al. 1241, carried out in the context of the chiral 
perturbation theory, which were used in extracting \V u d\ from the 7r e 3 decay, discussed earlier. Also, 
the so-called long-distance part of the electromagnetic radiative corrections, calculated by Ginsbe rg Ion, 
ago and used in the Leutwyler-Roos analysis EHI nas been recently checked (and corrected) l ^4 | 3' 
Finally, two 0(p 6 ) chiral perturbation the ory cal culations of the isospin-conserving contribution to the 
Kt3 form factor have also been undertaken EIE2I i n particular, it has been pointed out by Bijnens and 
Talaveral^l that the low energy constants (LEC's) which appear in order p 6 in the form factor /+(0) can 
be determined from measurements via the slope and the curvature of the scalar form- factor fo(q 2 ). 
In fact, there is some model-dependence also in the order p A parameters which impacts on \ V US \ : and one 
should firm up the existing phenomenological estimates by new measurements and/or calculations of the 
LEC's on the lattice. 

In addition to these theoretical developments, new experiments and/or analysis have been reported 
during this year by several groups. This includes a ne w, h igh statistics measurement of the K + — > tv e + v 
(K^ 3 ) branching ratio by the BNL experiment E865k™ which impacts on the determ ination of \V US \. 
New results in Kg^ decays have been reported by the KLOE collaboration at DA<i>NE OSES based on 
the measurements of the decays Kl — > Ttfiffj,, Kl —* 7r e^ e , and K$ — > Trei> e . In addition, semileptonic 
hyperon decays have been revisited by Cabibbo et al. ESI t determine the Cabibbo angle (or | Vu S \)- 
Finally, a determination of \V US \ has been undertaken from hadronic r-decays by Gamiz et al. "*-4 In 
this subsection, we summarize these results, some of which are new additions in this field since the 
CERN-CKM workshop. 

The four Kg^ decay widths for the decays K^ 3 , K® 3 , K~Z S , and K® 3 have been analyzed by 

Cirigliano et al. 1*^1 Normalizing the decay widths in terms of the quantity f? 71 (0), evalu- 
ated in the absenc e of th e electromagnetic corrections, the following master formula is used to ex- 
tract |K s |/f°^(0)ESE3 

o / r x 1 / 2 / i x 1 / 2 i 

\V us \f«* (0)= [tf nIn ( Xu0) ) [s^J i + ( jn [7(2)+( Jg M + A/ n (A t ,a)/2' (20) 

where the index n runs over the four modes, Af n = C 2 G 2 ? M^/192tt 3 , with C„ = (1, l/y/2) for (if , K + ). 

1771 

The various corrections and the compilation of the decay widths can be seen in the literature >=M This 
yields the following value: 

IK, I/+ (0) = 0.2115 ± 0.0015 . (21) 



The quantity fF * (0) has been studied in the context of the chiral perturbation theory. The result up 



i + 

to the next-to-next-to-leading order is known 

ff n ~(0) = l + h + U + O(p 6 ), (22) 

with f 2 = -0.023 and f 4 = -0.016 ± 0.008, yielding the Leutwyler-Roos value ff % ~ (0) = 0.961 ±0.008. 
This gives 

\Vus\Kia = 0.2201 ±0.0024. (23) 

Bijnens and TalaveraGH have included the isospin-conserving part of the 0(p 6 ) corrections in the determi- 
nation of ff^ (0), getting ff*~ (0) = 0.9760±0.0102, which, in turn, yields \V us \ Kt3 = 0.2175±0.0029. 
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However, as emphasized by these authors, this result should be treated as preliminary since the isospin- 
breaking 0(p 6 ) contributions are not yet included. Also, the effect of the curvature in the form factor on 
the experimental value remains to be evaluated. 

Recently, the E865 collaboration at Brookhaven ^31 has published a branching ratio for the decay 
K+ 3 ,y BR(K+ 3h] ) = (5.13±0.02(stat) ±0.09(sys) ±0.04(norm))%, which is about 2.3a higher than the 

current PDG valued for this quantity. The higher E865 branching ratio translates into a correspondingly 
higher value of the product \V us f+ n (0)|: 

\V US /f >O (0)| = 0.2239 ± 0.0022(rate) ± 0.0007(A t ) , (24) 

which on using the result from Cirigliano et al. I3BI for f£ Tr ° (0) yields 

|Ks|e865 = °- 2272 ± 0-0023(rate) ± 0.0007(A t ) ± 0.0018(/ + ) . (25) 

This differs from the older Kg 3 result 123|) by more than 2a. Interestingly, the E865 value of \V US \, 
together with the world average for \V u d\ given in (|19|) . leads to perfect agreement with the CKM- 
unitarity! Denoting the departure from unity in the first row of Vckm by Ai (defined earlier), the value 
obtained by the E865 group yields Ai = 0.0001 ± 0.0016^3. 

The other new addition to this subject is the measurements of \V US \ from the production of the 
^-meson at DA$NE and its decays into K + K~ and K^Kg pairs with the subsequent i^-decays. The 
KLOE collaboration at DA$NE will eventually measure \y us \ precisely (to an accuracy of better than 
1%) from all four channels of the K$ and Kl decays (involving the final states 7r liwl= = e, fi) as well as 
from the decays. Their preliminary results are available in conference reports I ' ' on the following 
three modes: K$ — > Ttei> e , Kl — > Tr/iv^ and Kl — > Ttei> e . Of these, the analysis of the K 5 decay mode 
is more advanced in terms of the systematics. Concentrating on this decay, its branching ratio has been 
measured by KLOE as BR{K S -> n-e + u e + c.c.) = (6.81 ± 0.12 ± 0.10) x 10~ 4 . The lifetime of the 
ifs-meson has been recently measured by the NA48 experiment t(K s ) = 0.89598(48)(51) x 10~ 10 s, 
allowing to have a new precise measurement of the ratio HR(Ks - * itev e )/T(Ks). Using the theoretical 
analysis of the K® 3 mode discussed earlier, this yields 

\V US I /f (0) = 0.2109 ± 0.0026 , (26) 

in excellent agreement with the value given in (|21|l . obtained from the earlier results on Ki 3 deca ys. The 
corresponding (preliminary) values from the two Kg 3 decay modes of the if^-meson are similar Iv^lvv l 
with |K s |/f°M0) = 0.2085 ± 0.0019 (for the nev e mode) and \V US |/f % ~ (0) = 0.2106 ± 0.0028 (for 
the niiis^ mode). However, as the systematic errors (in particular, for the Kl modes) have not yet been 
finalized, these numbers should not be averaged yet. Following the advice of the KLOE collaboration , 
we take the value of \V US \ obtained from the better studied K$ mode, as the preliminary value of this 
quantity from the current DA$NE measurements: 

IK, | D A*NE;/f s =0.2194 ±0.0030. (27) 

This is in comfortable agreement with the earlier determinations of this quantity given in (|23J) . 

While still on the subject of determining \ V US \ , there are two non-^3 estimates of this matrix element 
available in the literature, the first estimate is from the study of the semileptonic hyperon decays and 
the second is from the hadronic decays of the r-lepton. We discuss them in turn. 

The value listed in the PDG review for \V US \ from hyperon decays \V US \ = 0.2176 ± 0.0026 is similar 
in its precision as the one in (|23|l . obtained from the Kg 3 decays. However, based on the observation that 
the value obtained from hyperon decays is illustrative as it depends on the models to incorporate the 
S'?7(3)-symmetry-breaking corrections, and the theoretical dispersion (model-dependence) is significant, 
this value of |Ks| is not included in the world average of |Ks| by the PDG. This state of affairs was 
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where 



considered more or less as a theoretical fait accompli and no significant attempt was undertaken to reduce 
this model dependence. Recently, Cabibbo et al.l^have taken a somewhat different approach and have 
reported an analysis of the hyperon decays to extract \V US \. Their main assumption and results are 
summarized below. 

Denoting a ty pical hyperon decay by M\ — > A^e - !^, the following four decays are reanalyzed by 
Cabibbo et al.!2Sl (Af 1 ,JV 2 ) = (A,p) , (£~,n) , (H~, A) , (S°,£ + ). The matrix elements for these decays 
can be expressed as follows 

Gq 

M = —jL (AA 2 | J Q (0)|M) [uel a (l + 75H] , (28) 
v 2 

WMOm) = h(q 2 ) la + a a0 <f + 9l (q 2 ) la 75 + 9 -^fl <r a0 q? 75 , (29) 

Gs = G fV us (G fV ui i) for |A5| = 1(AS = 0) processes, and the contribution proportional to the electron 
mass has been dropped. The analysis by Cabibbo et al. focuses on the experimentally measured decay 
rates and the measured quantity <?i//i, which liberates them from estimating this ratio from theory. 
This is then used with the theoretical values of /1, f 2 , and 172, calculated in the SU(3)-symmetry limit 
to determine |V^ S |. Deviations from the SU(3)-symmetry limits of these quantities are expected to be of 
varying magnitude. Corrections to f\ are of second order, due to the Ademollo-Gatto theorem, but the 
weak magnetism f 2 is not protected by this theorem. Likewise, SU (3)-breaking effects invalidate the usual 
argument based on the absence of the second class currents and SU(3) symmetry, which yields g 2 = 0. No 
precise experimental information is available on g 2 . Expressing /2//1 in terms of the anomalous magnetic 
moments of the neutron and the proton, and appl ying the S't/(3)-symmetry to the ratio /2/Mjvi) where 
M/Vi is the mass of the parent hyperon, yields 

\V us \(A^pe-&) = 0.2224 ±0.0034, 

|K S I (E~ -► ne~u) = 0.2282 ± 0.0049 , 

(30) 

\V us \{zT ->Ae~D) = 0.2367 ±0.0099, 



giving an average value 



V us \ (S u -> S+e"P) = 0.209 ± 0.027 , 



|K S | H ypcron = 0.2250 ±0.0027. (31) 



While this analysis is internally consistent, namely that the values of \V US \ returned from the four decays 
are compatible with each other, and this observation is used by Cabibbo et al. ESI to argue that the data 
are compatible with the assumption that the residual SU(3)-breaking corrections are small, this feature is 
less transparent in model-dependent theoretical studies. It is difficult to quantify in a model-independent 
way the effects of 5J7(3)-breaking in fx and fa (as well as a non-zero value of (72)7 which are bound to 
renormalize the value of \V US \. Lattice calculations can clarify the theoretical issues involved, assuming 
that they will reach the required precision. Interestingly, the combined value of | V us \ from the hyperon 
decays in (|31|l together with the value of \V u d\ in l|19|) leads to perfect agreement with the CKM unitarity 
for the elements in the first row. 

Finally, we discuss the novel method advocated by Gamiz et al. E3 to determine \V US \ from the 
analysis of the hadronic decays of the r-lepton using the spectral function sum rules. In this method, 
\V U s\ and m s , the s-quark mass, are highly correlated and it is difficult to determine both. Since m s 
is known from other methods, one could fix its value in the current range, and optimise the analysis to 
determine | V us \ . This is what has been done by Gamiz et al. , which we briefly summarize below. 
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The starting point of this analysis is the moments Rf of the invariant mass distributions of the final 
state hadrons in the decay r — > v T + X: 

fl-.f 5 *^ «)'(^y*. (32) 

mi J \ mi J as 



Here R®° = R T , with R T defined as follows: 



R T = - „, = , vfi— = Rt,v + Rt.a + Rt.s 7 (33) 

where the vector (V), axial- vector (A) and scalar (S) contributions are indicated, with the scalar contri- 
bution coming essentially from the us branch of the decay r~ — > v T + us. 

The moments can be expressed in a form in which the dependence on IV^dj 2 and | Vu, s | 2 becomes 
explicit: 

R k T l =3(\V ud \ 2 + \V us \ 2 )A R h+**<(0) + £ (cos 2 c( 5^+sin 2 c C (I?) )j , (34) 

where the short-distance radiative correction Ar has been encountered earlier, 5 k " > is the perturbative 
dimension-0 contribution, and = (S^y + #y ^ )/2 stand for the average of the vector and axial 

vector contributions to the (kl) moments from dimension D > 2 operators in the operator product 
expansion of the two-point current correlation function governing r-decays. 

Theoretical analysis in the determination of \V US \ is carried out in terms of the SU(3)-breaking 
differences defined as: 

M -W^^ a \S(C--^>), (35) 

which do not involve the perturbative correction £ fc '(°) an d v anish in the SU(3) limit. Concentrating on 
the (0, 0) moment for the analysis, for which Gamiz et al. calculate SR T = 0.229 ± 0.030 for the r.h.s. 
of the above equation, using the experimental input EHI R T = 3.642 ± 0.012 and i? Tj s = 0.1625 ± 0.0066, 
and invoking CKM unitarity to express \V u d\ in terms of \V US \, yields 

iKslr-dccays = 0.2179 ± 0.0044(exp) ± 0.0009(th) = 0.2179 ± 0.0045 . (36) 

The first error is the experimental uncertainty due to the measured values of R T and R T .s, which is the 
dominant error at present but can be greatly reduced if the B-factory data on r decays is brought to bear 
on this problem, and the second error stems from the theoretical error in the calculation of 5R T , which 
is dominated by the assumed value for the s-quark mass: m s (2 GeV) = 105 ± 20 MeV, and should also 
decrease in future as the s-quark mass gets determined more precisely. While the current error on |V^ S | 
from T-decays is approximately a factor 2 larger at present than the corresponding error on this quantity 
from the analysis, potentially r-decays may provide a very competitive measurement of |K»s|- Of 
course, we also expect substantial progress on the front from the ongoing experiments. 

The present status of \V US \ is summarized in Fig.^ an d is based on the following five measurements: 
(i) From the old data, (ii) from the K^ 3 measurements by the BNL-E865 collaboration, (iii) from 
the KLOE data on K$ decays (still preliminary), (iv) from hyperon semileptonic decays, and (v) from 
r-decays. The current world average based on these measurements 

|Ks|wa = 0.2224 ±0.0017, (37) 

is also shown in this figure. We have added the statistical and systematic errors in quadrature. As not 
all the measurements are compatible with each other, we have used a scale factor of 1.3 in quoting the 
error. The resulting value of \V US \ is somewhat larger but compatible with the corresponding PDG value, 
IKJpdg = 0.2196 ±0.0026. 
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Figure 1. Present status of \V US \ measured in various experiments. Note that the KLOE entry is based on Kg decays. The 
band results from the current measurements of \V u d\ an d |V„f,|, and imposing the CKM unitarity constraint. 



2.3 Unitarity constraint for the first row in Vqkm 

In discussing the test of the CKM unitarity in the first row, we use the current world average of the 
matrix elements \V u d\ — 0.9739 ± 0.0005 to determine from the unitarity constraint a value for | V^ s |: 

I | un it = 0.2269 ±0.0024. (38) 

This is shown as a vertical band in Fig.^ The current world average |Kts|wA from direct measurements 
(|37|l differs from its value IKslu nit b y 1.5c. In this mismatch, \V u b\ plays no role, as its current value 
\V ub \ = (3.80l8;|| ± 0.45) x lO" 3 41 is too small. 

However, in averaging the value of \V U8 \, if one leaves out the entries from the BNL-E865 and Hyperon 
data, the former on the grounds of being at variance with the PDG value for BR(ii^), and the latter due 
to the neglect of the SU(3)-breaking corrections in some of the form factors, which can only be estimated 
in model-dependent ways, then the resulting world average goes down, yielding \V US \ = 0.2195 ± 0.0017. 
The error is almost the same as the one shown in (|37fl . as the remaining measurements are compatible 
with each other, and hence do not require a scale factor [S = 1). This value of \V US \, together with \V u d\ 
given above, yields Ai = (3.3 ± 1.3) x 10 -3 , corresponding to a 2.5<t violation of unitarity in the first 
row of the CKM matrix. Tentatively, we conclude that the deviation of Ai from zero is currently not 
established at a significant level. We look forward to the forthcoming results from the KLOE collaboration 
(as well as from NA48), on r-decays from the B-factories, and also improved measurements of Qa/ gv in 
polarized neutron beta-decays, which will yield more precise measurements of \V US \ and \V u d\, enabling 
us to undertake a definitive test of the unitarity involving the first row of the CKM matrix. 
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3 Current Status of \V c d\ and \V CS \ 

Concerning the determination of \ V c d\, nothing much has happened during the last decade! Current value 
of this matrix element is deduced from neutrino and antineutrino production of charm off valence d quarks 
in a nucleon, with the basic process being v^d — > /i~c, followed by the semileptonic charm quark decay 
c — > sfi + iSfi, and the charge conjugated processes involving an initial beam. Then, using the relation 

Q-(z^ -> - Cr(P M -> fl+fl') 3 +VMT/ 12 fort 

- f Z\ — = t;B[c^ V X )\V cd \ , 39 

one obtains \V cc t\ from the current average of the l.h.s., quoted by the PDG EI as (0.49 ± 0.05) x 10~ 2 , 
and B(c -> /J, + X), for which the PDG average is B(c -> (J.+X) = 0.099 ± 0.012, yielding 

| V cd \ =0.224 ±0.016. (40) 

Compared to \V US \, the precision on \V c d\ is not very impressive, with (5|U c d|/| V c d\ = 7%. 
Concerning \V CS \, three methods have been used in its determination: 

1. Semileptonic decays D — ► K£ + i>g, 

2. Decays of real W ± at LEP: W + -> cs(g) and -> cs{g), 

3. Measurement of the ratio T{W ± -» hadrons)/T(W ± -> ^ ± ^). 
We briefly discuss them in turn. 

|y c8 | from D -> 

This makes use of the following relation: 

r(D _ Ke ^ . ffi^fii _ SdM , l/+(fl)l . (1 + &) , (41) 

TD 1927T' 3 

where $ is the phase space factor, /+(0) is the dominant form factor in the decay evaluated at 
q 2 = (requiring extrapolation of data to q 2 = 0), and Sr arises from the g 2 -dependence of this form 



factor. Using ™* /+(0) = 0.7 ± 0.1, coming from the early epoch of the QCD sum rules, the earlier 
version of the PDG CKM review E2l quotes a value \V CS \ = 1.04 ± 0.16. The decays D -> Kiv t and 
D — > K*tvi have also received quite a lot of attention by the lattice groups in the past. There is almost 
a decade old result from the UKQCD collaboration 1441 /+(0) = 0.67^0 08' an< ^ a relatively recent result 
by the Rome Lattice group ESI, /+(0) = 0.77 ± 0.041 q'q 1 . In fact, the lattice technology has advanced 
to a point where precision calculations of the relevant form factors in D — > K and D — > K* can be 
undertaken. On the experimental side, one already has a measurement of the form factor in D — > K, 
and ratios of the form factors in D — > K*lvg have now b een measured with good accuracy, most recently 
by the FOCUS photoproduction experiment at Fermilab'-^. However, in the current version of the PDG 
review"^, the determination of y/ cs \ from D — > K transitions has been dropped. We hope that in fu- 
ture, with improved theory and experiments, this value judgement on the part of the PDG will be revised. 

\V CS \ from W + — ► cs(g) and W~ — > cs(g): 



This method involves the process e + e — > W + W , well measured at CERN, and subsequent 
charmed-tagged W^-decays. The ratio 



Acs) _ 



T(W+ -► cs) /T(W+ -> hadrons) , (42) 
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then allows to determine |T4 S |. The weighted average of the ALEPH I^ZI and DELPHI ESI measurements 
yields^ 

\V CS \ = 0.97 ± 0.09 (stat.) ± 0.07 (syst.) . (43) 
The current precision on \V CS \ from direct W -measurements is 5| V^ S |/|V^ S | = 11%. 

Measurement of the ratio T(W ± — > hadrons)/T(W ± — ► £ + vg): 

A tighter determination of |T4 S | follows from the ratio of the hadronic W decays to leptonic decays, 
which has been measured at LEP. Using the relation 

yields Ei 

\Vij\ 2 = 2 - 039 ± °- 025 (B(W -> ±0.001 (a s ), (45) 

and gives (on using the known values of the other matrix elements) 

\Vc\ = 0.996 ±0.013. (46) 

The measurement ()45f) provides a quantitative test of the CKM unitarity involving the first two rows of the 
CKM matrix. This amounts to a violation of unitarity by 1.6cr, and hence is statistically not significant. 

The matrix elements \V c d\ and \V CS \ will be measured very precisely in the decays D — > Kive and 
D — > n£v£ by the CLEO-C and BES-III experiments, with anticipated integrated luminosity of 3 fb~ 
and 30 fb _1 at ^(3770), respectively. These experiments will also allow, for the first time, a complete 
set of measurements in D — > (K, K*)£vg and D — > (n, p)lvi of the magnitude and slopes of the form 
factors to a few per cent level. From a theoretical point of view, this is an area where the Lattice-QCD 
techniques can be reliably app lied to enable a very precise determinations of the matrix elements \V CS \ 
and \V c d\. Typical projections ^21 at the CLEO-C are: <J|V^ S |/|T4 S | = 1-6% and a similar precision on 

s\v cd \/\v cd \. 

4 Present Status of |T4&| and \V u b\ 

The matrix elements \V c b\ and \V u b\ play a central role in the quantitative tests of the CKM theory in 
current experiments. In particular, these matrix elements enter in the following unitarity relation 

v ud v: b + v cd v; b + v td v; b = o . (47) 

This is a triangle relation in the complex plane (i.e. p-fj space). The three angles of this triangle are 
defined as: 

( v t * b v td \ a ( v: b v cd \ ( v: b v ud \ 



The BELLE convention for these phases is: <p2 — a, <fii — (3 and 03 = 7. In the Wolfenstein parametriza- 
tion given above, the matrix elements V ud , V cd , V cb and V tb entering in the above relations are real, 
to 0(A 3 ). Hence, the angles (3 and 7 have a simple interpretation: They are the phases of the matrix 
elements Vt d and V ub , respectively: 



Vtd = \Vtd\e~ lfi , V ub = \V ub \e~^ ; (49) 
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and the phase a defined by the triangle relation: a — ir — (3 — 7. The unitarity relation 147|) can be 
written as 



JW 7 + Rte~ l(3 = 1 



(50) 



where 



^6 




vV + n 2 = ( 1 



2 



) 



A V cb 



1 V ub 



(51) 



\VttVtd\ 

\v; b v cd \ 



A V cb 



1 V td 



The unitarity triangle with unit base, and the other two sides given by R b and R t , and the apex defined 
by the coordinates (p, fj) is shown in Fig. [21 Quantitative tests of the CKM unitarity, being carried 
out at the B factories, consists of determining the sides of this triangle through the measurements of 
\V ub \, \Vcb\ and |Vt<j| as precisely as possible, which allows to determine indirectly the three inner angles 
a, 0, 7, and confronting this information with the direct measurements of the three angles (a, /?, 7) (or 
4>U 4>2-> fa) through the CP-violating asymmetries. We will return to a quantitative discussion of these 
tests in Section 6. 



Figure 2. The unitarity triangle with unit base in the p - r) plane. The angles a, /3 and 7 are defined in 1481 and the two 
sides ijj, and Rt are defined in 1511 . 

Current status of the matrix elements \V cb \ and \V ub \ has been discussed in the literature in great 
detail. These include proceedings of the research workshops and con ferences on flavour phys ics held 
rccc in particular, the experimental reviews by Gibbons^, ThorndikejUj and Calvi^, and 

theoretical developments reviewed by Luke^l Ligeti^^, LellouchE3 and Uraltsev^l Updated data are 
available on the webcites of the working groups established to perform the averages of the experimental 
results in flavour physics ESI and the CKM matrix elements ESI. We shall make extensive use of these 
resources, focusing on some of the principal results achieved, and discuss their theoretical underpinnings. 

4-1 Present status of \V cb \ 

Measurements of \V cb \ are based essentially on the semileptonic decays of the 6-quark b — ► clvg. In 
the experiments, one measures hadrons, and hence the inclusive hadronic states B — > X c li>t and some 
selected exclusive states, such as B — ► (D, D*)£vg, is as close as one gets to the underlying partonic 
weak transition. In the interpretation of data, QCD is intimately involved. In fact, quantitative studies 




(p, v) 




13 



of heavy mesons, in particular B-mesons, have led to novel applications of QCD, of which HQET 
in its various incarnations is at the forefront. Scmileptonic £?-decays have also received a great deal of 
theoretical atte ntio n in methods which involve non-perturbative techniques, foremost among them are the 
QCD sum rules 1^ and Lattice QCD We shall restrict the theoretical discussion to these frameworks. 



4-2 Determination of \V c b\ from inclusive decays B — > X c iv(, 

The theoretical framework to study inclusive decays is based on the operator product expansion (OPE) , 
which allows to calculate the decay rates in terms of a perturbation series in a s and power corrections 
in Aqcd/w& (and Aqcd/w c ). This tacitly assumes quark-hadron duality, which is supposed to hold 
for inclusive decays and also for partial decay rates and distributions, if summed over sufficiently large 
intervals (3> Aqcd), and weighted distributions (moments). Deviations from this duality are, however, 
hard to quantify, and they will be the limiting factor in theoretical precision ultimately. The first term of 
this QCD corrected series is the parton model result for the decay b — > cinug. Leading Q(a s ) correct ions 
were obtained some time ago for inclusive decay rates EH and lepton energy distribution ESESE3 I n 
the meanwhile, the QCD perturbative corrections to the decay rates are known up to order aj.(3o 
where (3q = 11 — 2n//3, with rij being the number of active quarks. This term usually dominates the 
0(a 2 s ) corrections, though there exists at least one counter example, namely the inclusive decay width 
T(B — > X s j), where the contribution to the width in a 2 /3o is small ESI. The result in the MS scheme is EH1 



T(b^ceP e (+g)) = \V cb \"^4o.52 

19z7T 



1_ 1.67(^1) - 15.1 



(52) 



where the numerical coefficients correspond to the choice m c /m& = 0.3. 

The leptonic and hadronic distributions in B decays are now calculated using techniques based on the 
OPE. While the distributions themselves are not calculable from first principles and invariably involve 
models, called the shape functions, inclusive decay rates, partially integrated spectra, and moments are 
calculable in the OPE approach in terms of the matrix elements of higher (than four) dimension operators. 

The book-keeping of the power corrections is as follows. The leading A/ mi, correction to the decay 
rates vanishes in the heavy quark li mit ^4 and t he 0(A 2 /m\) effects can be parametrized in terms of 
two parameters Ai and A2, defined as I 71 | 72 | 73 I 74 | 



Ai = -^—(B(v)\b v (iD) 2 b v \B(v)) , A 2 = J—(B(v)\b v ^-a^G^b v \B(v)) , (53) 



where b v denotes the b quark field in HQET, with D M and being the covariant derivative and the QCD 
field strength tensor, respectively. While Ai is not known precisely, having a value typically in the range 
Ai = [-0.1,-0.5] GeV 2 , but A 2 is known from the B* - B mass difference to be A2(to&) — 0.12 GeV . 
Corrections of order a s jm\ are still not available, but 0(1/ mf) corrections to the decay widths have 
been calculated They are expressed in terms of six additional non-perturbative parameters, pi, P2, 
and 71, i = 1, 4. There are two constraints on these six parameters, which reduces the number of free 
parameters relevant for B-decays in this order to four. So, including the 0(l/m^) terms, there are in 
all six non-perturbative parameters which have to be determined from experimental analysis. They can 
be determined, or at least bounded, from the already measured lepton- and hadron-energy moments in 
B — > X c li>i and photon-energy moments in the inclusive decay B — > X s j. 

Theoretical results for the inclusive rate and moments depend on the scheme for defining the b- 
quark mass. They influence the decay rates more, in particular the branching ratio B(B — > X s "f), 
where the sch eme- dependence of m,b and m c is currently the largest theoretical uncertainty, typically of 
order 10%E§ but the moments are less sensitive. We shall discuss here the results in the so-called T(15 l ) 
schemeEZEHl 

{m\ s denotes the &-quark mass in this scheme and my is the T-meson mass), as this scheme 
is en vogue in the analyses of the moments by experimental groups. Taking into account the perturbative 
corrections to order ce^Po and power corrections to order l/mj), the result for the semileptonic decay 
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width T(B — > X c lvi) in T(IS') mass scheme reads as follows 



E2 



G F \V cb \ 2 /TOTN 5 

192tt 3 V 2 



0.534 - 0.232 A - 0.023 A 2 - 0.11 Ai - 0.15 A 2 - 0.02 AiA + (54) 



+ 0.05 A 2 A - 0.02 pi + 0.03 p 2 - 0.05 T x + 0.01 T 2 - 0.07 % ~ 0.03 T 4 - 0.051 e - 0.016 e| LM + 0.016 eA 

where e and £blm are parameters in the perturbative part, entering through the a 2 /3o term EE1 the 
parameter denoted by A is called A = rag — in the HQET jargon and is defined above in terms of the 
T(15) mass by A = my/2 — m\ . The corresponding expressions for the decay widths in other schemes 
for the 6-quark mass can be seen in the paper by Bauer et al. 

Determinations of \V c b\ from this method are based on the analysis of the following measures. For the 
lepton energy spectrum, partial rates and moments are defined by cuts on the lepton energy (E{ > Eq,Ei): 

d /El 77i \ / Ei dEp dEf, J g o Ef, ^ dEi 

KoiEo,^)-- SF — , u n (E )-— dr , (55) 

Je„ dEl dbj t Jeo dE- e dh e 

where dT/ dEi is the charged lept on spectrum in the B rest frame. The moments R n are known to 
order a 2 /?o and Aq CD /to 3 ^51. For the hadronic moments, the quantities analyzed are the mean 
hadron invariant mass and its variance, both with lepton energy cuts Eq: Si(Eq) = (n^x ~ m2 n)\E,>E n 
and S 2 {E ) = {{m x — ( , m x )) 2 }\ Ee>Eo , where m D — (m D + 3m_o*)/4. S n are known to order a 2 /3 
and Aq CD /m 3 Lu4 For the decay B — + X s j, the mean photon energy and variance with a photon 
energy cut, £7 7 > Eq, calculated in the B rest frame, have been used: Ti(Eq) = (E~ / )\e^>e and 
T 2 (E ) = ((E 7 - (£ 7 )) 2 >| B .,>£ . Also, T 1>2 are known to order af/3 ^and A| CD /m 3 ESl. 

This theoretical framework has been used by the CLEO 84 , BABAR — and DELPHI EH1 collabora- 
tions, and their results for \V c b\ using the moment analysis are as follows: 

\V cb \ = (40.8 ± 0.6 ± 0.9) x 10~ 3 [CLEO] , 

\V cb \ = (42.1 ± 1.0 ± 0.7) x 10~ 3 [BABAR] , (56) 

\V cb \ = (42.4 ± 0.6 ± 0.9) x 10 -3 [DELPHI] . 

The moments (and hence the values for \ V c b\) are strongly correlated with m,b, and are scheme-dependent. 
This aspect should not be missed in comparing or combining various determinations of \V cb \- They can 
be averaged to get |V^,| from inclusive decays 

iKblind = (42.1 ± 0.7 cxp ± 0.9 thco ) x 10~ 3 , (57) 

where we have kept the theoretical error as 0.9, which is an underestimate as no account is taken of the 
duality error, probably not negligible at this prec ision. A value very similar to the ab ove results was 
obtained by Bauer et aJ,. E5I using the CLEO dataEZEBEHl the earlier BABAR dataEHl and data from 
the DELPHI Collaboration ED 

\V cb \inol = (40.8 ± 0.9) x 10~ 3 [Bauer et al.} . (58) 

A mismatch between the earlier BABAR measurements E01 f the hadron invariant mass spectrum 
presented as a function of the lepton energy cut {Eq) and the corresponding OPE-based theoretical 
analysis of Bauer et al. is now largely gone. The updated BABAR ESI anc j CLEO E3I data are in 
agreement with each other and with the OPE-based theory. This is depicted in Fig. E] showing the 
hadron moment (M x — Mf>) vs. lepton energy cut. Theory bands taking into account the variati ons in 
the input parameters are also given and the details of the analysis can be seen in the CLEO paper E2I 

Mutual consistency of the experiments in terms of \V c b\ and m I is show n in Fig. Th e contours 
represent the best fits (A% 2 = 1) for the hadron moments (from the BABAR ESI CLEoEBl anc i DEL- 
PHI ESI data) and lepton moments (from the CLEO El and DELPHI ESI data). One observes from these 
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Figure 3. Comparison of BABAR '03 CLEO '03 EH and CLEO '01 measurements of the hadron invariant mass 
moment (M^ — M^) vs.lepton energy cut. The bands show the parametric uncertainties from theoryEHl (Figure taken 
from the CLEO paperE3.) 



correlations that there is still some residual difference between the best fit contours resulting from the 
analysis of the lepton- and hadron-energy moments. The current mismatch remains to be clarified in 
future experimental and theoretical analyses. Once the experimental issues are resolved, the {\V c b\ - 
m\ s ) correlation from the hadron and lepton moments can be used to quantify the quark-hadron duality 
violation. 

4-3 Determination of \V c b\ from exclusive decays B — > {D,D*)£vi 

In exclusive decays, B — ► (D, D*)£i/g, one needs to know the hadronic matrix elements of the charged 
weak current, {D\J^\B) and (D*\J^\B). The former involves two form factors, called F + (q 2 ) and F (q 2 ), 
and for the transition to the _D*-meson, one has four such form factors, called in the literature by the 
symbols V(q 2 ), Ai{q 2 ), i = 1, 2, 3. If these form factors can be measured over a large-enough range of q 2 
and the same can be obtained from a first principle calculation, such as lattice QCD, then exclusive decays 
would provide the best determination of \V c b\- In the absence of a first principle calculation of these form 
factors, HQET provides a big help in that the heavy quark symmetries in HQET allow to reduce the 
number of independent form factors from six in the decays at hand to just one, called the Isgur-Wise 
(IW) function l^sl ^(w), where u — v.v', with v and v' being the four-velocities of the B and D(D*) 
meson, respectively. Moreover, HQET provides a normalization of the IW function at the sym metry 
point, u = 1. A lot of attention has been paid to the decay B — > D*iv^ due to Luke's theorem ^3 which 
states that symmetry-breaking corrections to Tiu = 1) are of second order, a situation very much akin 
to the Ademollo-Gatto theorem for the form factor /+(0) discussed earlier. 
The differential decay rate for B — * D*£vi can be written as 

£ = (c 2 - 1) 1/2 ml, (m B m D « ) 2 Q{ui) \V cb \ 2 \T(uj)\ 2 , (59) 
where G{lS) is a phase space factor with <?(!) = 1. Theoretical issues are then confined to a precise deter- 
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Figure 4. Constraints on the HQ ET pa ram eters obtained from hadronic moments as measured by BABAR ES| hadronic 
moments combined (BABAR EST CLEoEHl and DELPHI EHJ) and the combined lepton energy moments (CLEO Inland 
DELPHI EHI) in the (\V cb \,ml s ) plane. (Figure taken from the BABAR paper ESI.) 



mination of the second order corrections to T{io = 1), the slope of this function, p 2 , and its curvature c, 

= [1 + P 2 (oj - 1) + c(oj - l) 2 + ...] . (60) 

In terms of the perturbative (QED and QCD) and the non-perturbative (leading 5\/ m 2 and subleading 
<5i/m 3 ) corrections, the normalization of the Isgur-Wise function can be expressed as follows: 



= rj [l + Si/ m a + S 1/r , 



(61) 



wher e 77 is t he perturbative renormalization of the Isgur-Wise function, known to two loops, rj — 0.960 ± 



The formalism for calculating Si/ m 2 corrections in HQET has been developed by Falk and 



0.007 

Neubert^Hland Mannel^. The various non-perturbative parameters entering in 8i/ m a and the slope p 2 
have been stud ied in the context of quark models^^, sum rules ^ ^ ^^ j ^ ^^^^ l and quenched Lattice 
q CD | 106 | 10 Z| The default value used in the BABAR Physics BookGHEl jr(i) = Q.91± 0. 04 h as recently 
been confirmed in the quenched lattice QCD calculations, including also <5i/ m 3 corrections! I The slope 
p 2 is obtained by a simultaneous fit of the data for ^"(l)]^] and p 2 , and experiments use a form for T{u>) 
given by Caprini et al. IliliJI The resulting values of j V^b| and the ^ r (l)|T4h| - p 2 correlation from a 
large number of measurements from the LEP experiments, CLEO, BELLE and BABAR are summarized 
in Fig. EI They lead to the following world averages 



T(l)\V cb \ = (36.7 ±0.8) x 10" 



p 2 = 1.44 ±0.14 [ X 2 = 30.3/14] , 



(62) 
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ALEPH 

33.6 ± 2.1 ± 1.6 
OPAL (partial reco) 

38.4 ± 1.2 ± 2.4 
OPAL (excl) 
39.1 ± 1.6 ± 1.8 
DELPHI (partial reco) 
36.8+ 1.4 ± 2.5 
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36.7 ± 1.9 ± 1.9 
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43.6 ± 1.3 ± 1.8 
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Average 

36.7 + 0.8 
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Figure 5. Present status of .F^lJIVcj,! (left frame) and the J 71 (1)1 ^cE> I ~ P 2 correlation (right frame) from B — * D'lvi dec ays. 
Note that F(l) is called T{\) in the text. (Figures taken from the Heavy Flavor Averaging Group [HFAG LP2003] ESI.) 



which with the value .F(l) = 0.91 ± 0.04 leads to 

\V cb \ B ^ D *eu e = (40.2 ± 0.9 cxp ± 1.8 theo ) x 10~ 3 . 



(63) 



This value is in very good accord with the determinations of \V c b\ from the inclusive decays B — > X c lvg 
given in (|56(1 for the T(lS')-scheme. While this consistency is striking, the agreement among the various 
experiments in the ,F(1)|14&| - p 2 correlation from the decay B — > D*lvi is less so, having a rather high 
X 2 j X 2 /d.o.f. = 2.16. A robust average for \V c b\ from both the inclusive and exclusive decays is not yet 
available from the Heavy Flavor Averaging Group HFAG 1^1 It is also not clear to me how to do this as 
\V c b\ from the inclusive measurements is obtained in the T(15) scheme, as the quark masses have been 
defined in this scheme in the analysis of data, whereas \V c b\ from the exclusive decays does not have this 
dependence. An average can be given by weighting the two measurements with their experimental errors 
only, leaving the theoretical errors as they are, or one could add them in quadrature assuming they are 
independent. This gives 

\Vc\ = (41.2 ± 0.8 cxp ± 2.0 t hco) x lO- 3 = (41.2 ± 2.1) x 10~ 3 , (64) 

yielding a precision <5 1 V^j, | / 1 V^t, | ~ 5%. In view of the still open issues (such as duality-related error, quark 
mass scheme-dependence in the inclusive decay, large x 2 /d.o.f. in B — > D*tvi decay, which makes the 
various experiments compatible with each other only at the expense of an increased error, etc.), a more 
precise value of |V^(,|, in my opinion, is not admissible at present. However, despite all the caveats, the 
achieved precision in \V c b\ is indeed remarkable. 



4-4 Present status of \V u b\ 

Essentially, there are two methods to measure \V u b\- The first one analyses the inclusive decays B — > 
X u ivi 1 for which the branching ratio lies about a factor 60 below the dominant process B — > XJLvi. This 
circumstance makes it mandatory to apply harsh cuts to tag well the b — > u events, invariably bringing in 
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its wake theoretical problems involving enhanced non-perturbative and perturbative effects. The second 
method involves the exclusive decays, such as B — > (n, p, uj)£u£, which require good knowledge of the 
form factors, not yet completely under theoretical control. We briefly summarize the present status for 
both the inclusive and exclusive determinations of |Vu&|. 



4-5 \V u b\ from inclusive measurements 

The theoretical framework to study the inclusive decays B — > X u lvg is also based on the ope rator product 
expansion. Up to 0(\/m 2 ) order, the result for T(B — > XJtvg) can be expressed as foilows^EB 

T(B -> X M = G ^ V f ( 1 - 2.41 £5sW - 21.3 (^Y + + , (65) 

1927T' 3 \ 7T V 7T / 2mf i / 

of which the first three terms are coming from the perturbative-QCD improved parton decay b — > utvg. 
The largest uncertainty in the decay rate is due to the 6-quark mass, mj, which is also scheme depe nde nt, 
as already discussed. In the T(lS')-scheme, the result for \V U b\ is numerically expressed as follows^!! 

, , , N o fB(B~*X u £v) 1.6ps\ 1/2 , N 

|K 6 | = (3.04±0.06 (pcrt )±0.08 ( . m6 ))xl0- 3 ^ v ^ ' J j , (66) 

where the first error has a perturbative origin and the second is from Amj, in the T(15)-scheme, ml = 
(4.73 ± 0.05) GeV. If this rate can be measured without significant cuts, then \V u b\ can be measured 
with an accuracy of 0(5%). We shall take this as an ideal case and discuss now the realistic cases when 
kinematic cuts are imposed on some of the variables to measure the b — > u semileptonic decays. 

As already mentioned, the dominant background is from the decays B — » X c ivg. Noting that the 
lowest mass hadronic state in B — > X c lvg is the D-meson, thus its mass mo is used to define the cut 
region to suppress b — > c transitions. So, the kinematic cut is either (i) on the upper end of the lepton 
energy spectrum, with Eg > {m 2 B — mp)/2mg, or (ii) on the momentum transfer squared of the tvg pair, 
q 2 > (rns — iri£)) 2 , or (iii) on the hadron invariant mass, mx < mp, or (iv) an optimized combination 
of some or all of them. These cuts reduce the experimental rates for B — * X u £ug, a handicap which will 
be overcome at B-factories with O(10 8 ) BB mesons already at hand. However, the cuts also make the 
theoretical rates less rapidly convergent in terms of the perturbation series in Aqcd/w{, and a s (mb). The 
other disadvantage is that the theoretical rate with a cut depends sen sitive ly (except for a cut on q 2 ) on 
the details of the S-meson wave function, or the shape function /(fc + )^ni. Here k + is the + component 
(using light cone variables) of a residual momentum of order Aqcd, entering through the relation 
= m^v^ + k^, where is the momentum of the 6-quark in the B meson, and i> M is the four- velocity 
of the quark. This can be seen as follows. With either of the two cuts, Eg > (m 2 B — m 2 D ) /2ms, or for the 
small hadronic invariant mass region, mx < mo, we have 

E x ~m b ; m\ — (m b v - q) 2 + 2E x k + + (67) 

bringing in the dependence on k + . This dependence is rather mild using the cuts on q 2 . The effects 
of the kinematic cuts on the decay distributions and rates have been studied at great length in the 
literature. I n fact, this enterprise has led to a flouris hing industry - the (kinematic) cutting technology 
using HQETm ^ 1 ^ 114 ! 11 ^ 116 ! 11 ^ 11 ^ 11 ^ 12 ^ 121 ! Some applications are discussed here. 

In the leading order in Aqcd/w;,, there is a universal function which governs the shape of the charged 
lepton energy spectrum, the hadronic in varia nt mass spectrum in B — » X u ivg and the photon energy 
spectrum in B — > X s j, defined as follows E2 

f(k+) = -^—(B(v)\b v 6(k+ + iD.n)b v \B(v)), (68) 
2m B 

where n is a light-like vector satisfying n.v = 1 and n 2 = 0. The physical spectra are obtained by con- 
voluting the universal shape function with the perturbative-QCD expressions. Ignoring the perturbative 
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and subleading power corrections, a measurement of the photon energy spectrum is a measurement of 
the shape function f(k + ): 

""' —{B^X u lv l )=[du6{m b -2Ei-u)f{u) + ..., (69) 



2r (°) dEi 

2fi-{B -» x sl) = J dujd{m b - 2£ 7 - u) f(u) + ... = f(m b - 2E- ( ) + . . 

where the normalization constants for the semileptonic and radiative B-meson decays are: 

(0) G 2 F \V ub \ 2 ml _ 2 
^ = 192^ = C *\V*\ > 



if = Uf ^ZT mb \cf\ 2 - \v tb v; s \\ (70) 



Kb 



32tt 4 

and Cf is the effective Wilson coefficient governing the decay -B — > A s 7. Thus, combining the data 
on _B — > X s 7 and the lepton energy spectrum from £> — > X u £ve, one can determine in the SM the 
following ratio: 

,2 =3«\Cf\ 2 ^ + 0(a s )+o(^l), (71) 
7r T S {E C ) \ m b J 

where T U {E C ) and T S (E C ) are the cut-off dependent decay width in T(B — > X u £vf, Eg > i? c ) and the 
cut-off dependent first moment of the photon energy spectrum, respectively 

V 1 Je c dE e ' 

2 r mB / 2 dV 
T S (E C )^— dE y (E^ — E c ) —2- . (72) 

m b J Ec dE 1 

It should be stressed that the ratio H71J1 holds not only in the SM, but also in models where the flavour 
changing (FC) transition b — > s is enacted solely in terms of Vckm, such as the minimal flavour violating 
supersymmctric models. An example where this relation does not hold is a general supersymmetric model 
in which the couplings diSjg, involving a down-type quark, a squark and gluino, are not diagonal in the 
flavour (ij) space. In that case, the decay width for B — > X s j does not factorize in |V t ;,V^* | 2 and depends 
on additional FC parameters. 

The relation (|71|) has been put to good use, in the context of the SM, by the CLEO collaboration H2i3 
through the measurement of the photon-energy spectrum in B — > A^l^land the lepton energy spectrum 
in B — > X u £i/£, with 2.2 GeV < E% < 2.6 GeV, as the photon energy spectrum has been well measured 
in the overlapping range for _E 7 , yielding 



\V ub \ = (4.08 ±0.34 ±0.44 ±0.16 ±0.24) 



xl0~ 3 , (73) 



where the first two uncertainties are of experimental origin and the last two are theoretical, of which 
<5|T4f>| = ±0.24 x 10 -3 is an assumed uncertainty in the relation H71JI. Combining all the errors in 
quadrature leads to \V ub \ = (4.08 ± 0.63) x 10 -3 , which is a ±15% measurement of this matrix element. 

This method of determining \V ub \ pioneered by the CLEO collaboration has received a lot of theo- 
retical attention lately. In particular, the subleading-twist contributions to the lepton and photon energy 
spectr a in the decays B — > X u tvi and B — > X s j, respectively, have been calculated in a number of 
papers H 24 | 125J126 | 127| providing estimates of the subleading correction in H71JI indicated as 0( A ^ p ). 
An important feature emerging from these studies is that the various spectra are no longer governed by 
a universal shape function f(u>), a feature which is valid only in the leading twist. In subleading twist, 
HQET shows its rich underlying structure leading to a number of additional subleading shape functions, 
which are no longer universal. The operators Oi needed to calculate the subleading twist contributions 
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and the corresponding matrix elements (shape func tions) of these operators (B\Oi\B) can be found, for 
example, in the papers by Bauer, Luke and Mannel ^^^^ . The photon energy spectrum in B — » X s "f 
can now be expressed in terms of three structure functions F(u>), h\{ui) and ffaC^) as follows 



m b dT _ j 



If \ . / A 2 



(4£L - m b )F(m b - 2£L) + Ui(m b - 2£ 7 ) + i7 2 (m b - 2E 1 ) + 

771 f, 1 ' 



QCD 



(74) 

where C 7 has been defined earlier. Here, F(tu) contains both the leading and sub-leading parts with 
F(ui) = f(ui) + 0(Aqcd/™&), and hi(u>) and -ff 2 (w) are the subleading shape functions. 

The corresponding lepton energy spectrum in the decay B — > X u lv( now has the following form^^ 



m b dT 2 
\Vub\ 



2C*si dE e 



dLu6(m b -2E t - u) ( F(u) I 1 - — ) - —hi(w) + —H 2 (lo)] +0 ( 



(75) 

where C s \ has also been defined earlier. It is obvious that the measurement of either the photon energy 
spectrum or the lepton energy spectrum does not allow to determine all three shape functions. Hence, 
they will have to be modeled . Th ese subleading corrections modify the relation Q71f) used in extracting 
\Vuh\, which can be written a 



V ub 



3 « |C e ff|2 rV^)y /2 (1 + (5(i?c)r (76) 



v tb v t * V ' 1 r s (£ c ), 

Again, S(E C ) can only be estimated in a model-dependent way. Typical estimates are S(E C — 2.2 GeV) ~ 
0.15, with 8{E C ) decreasing as the lepton-energy cut E c decreases, estimated as 0(10%) for E c — 2.0 GeV. 
One should use the order of magnitude estimate of the subleading twist contribution S(E C ) to set the 
theoretical uncertainty on \V ub \ from this method, which typically is 15%. 

These uncertainties can be reduced if one considers more complicated kinematic cuts, such as a 
simultaneous cut on mx and q 2 | H^ | 55 | wnose effect has been studied using a model for the leading- 
twist shape function f{k + ). The sensitivity of the partial decay width T(q 2 > q 2 ut ,mx < m cut ) on 
f(k+) is found to be small, and this is likely to hold also if the subleading shape functions are included. 
This method of determining \V ub \ has been applied by BELLE using two techniques. The first uses the 
decays B — > D^*'ivt as a tag, and the other uses the neutrino reconstruction technique, as in exclusive 
semileptonic decays, combined with a sorting algorithm (called "annealing") to separate the event in a 
tag and a b -> ulv g side. The method based on £>M-tagging yieldJi 2 ^ | V ub \ = (5.0 ± 0.64 ± 0.53) x 10~ 3 . 
The result using the annealing method with the cuts mx < 1.7 GeV, q 2 > 8 GeV 2 is ^^ l 

\V ub \ = (4.66 ± 0.28 ± 0.35 ± 0.17 ± 0.08 ± 0.58) x 10~ 3 , (77) 

where the errors are statistical, detector systematics, modeling b — > c, modeling b — > u, and theoretical, 
respectively. Combining all the errors yields \V ub \ — (4.66 ± 0.76) x 10~ 3 . 

A similar analysis by the BABAR collaboration, in which one of the two JS-mesons is constructed 
through the hadronic decays B — > D^h, and the inclusi ve se mileptonic decay of the other B-meson is 
measured with the cuts Ei > 1 GeV and mx < 1-55 GeV U^Ol yields 

\V ub \ = (4.52 ±0.31 ±0.27 ±0.40 ±0.25) x 10~ 3 , (78) 

where the errors are statistical, systematics, due to extrapolations to the full phase space, and from the 
HQET parameters, respectively. 

Finally, the effects of the so-called weak annihilation (WA)E3 which are formally of 0(AQ CD /m b ) 
but are enhanced by the phase space factor 167r 2 (compared to that of b — > uiui), introduce an additional 
theoretical uncertainty! I They stem from the dimension-6 four-quark operators in the OPE, 

Ov-A = (b v ^P L u)(u^P L b v ) , Os-p = (b v P L u){uP R b v ) , (79) 
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where Pl,r = (1 ± 7s)/ 2. In the lepton energy spectrum from B — ► X u lvt, they enter as delta functions 
near the end-point I f 

d ^ = -^f^f B rn B{Bl -B 2 )S(l-y), (80) 
ay lZir 

where y = 2Ei/irib and Jb — 200 MeV is the B meson decay constant; B\ and B2 parameterize the 
matrix elements of the operators Ov-A and Os-p, respectively: 

1 -{B\O v . A \B) = ^-B 1 , J-(B\Os- P \B) = ^B 2 . (81) 



2m B 8 2m B ' 1 1 ' 8 

In the vacuum insertion approximation, i.e., assuming factorization, their effect in the spectrum vanishes, 
as in this approximation, B\ = B2 = 1(0) for the charged (neutral) B mesons. Hence, they are generated 
by non-factorizing contributions and are not yet quantified. The se matrix elements are also encountered 
in calculating the difference s in the B^ and B° lifetimes E23 anc [ we refer to a recent discussion in 
the context of lattice QCD E3I However, comparing the extraction of \V u b\ from B ± — > X u lvi and 
B° — ► X u £i/g near the end-point of the lepton energy spectrum, one can determine the si ze of the WA 
effects. For the B^ decays, they can also be estimated from a related process B^ — > ^ 1/^7 ESI. 

The current results on \V u b\ from various inclusive measurements by the LEP, CLEO, BABAR and 
BELLE experiments are summarized by HFAgEDI No averaging for \V u b\ has been undertaken so far by 
this working group. However, the results in (JZHJ, O and (EHJ from the CLEO, BELLE and BABAR 
collaborations, respectively, have been averaged by Muheim.1 **" lto get \V u b\ from the T(45) data, yielding 

iKblinci = (4.32 ± 0.57) x 10" 3 , (82) 

in agreement with the LEP averageEZI \y ub \ = (4.09 ± 0.70) x 10~ 3 . 

4-6 \V u b\ from exclusive measurements 

First measure ments of the exclusive decays B — ► nlvi and B — » plvi were reported by the CLEO c ollab- 
oration in 1996^^1 Improved measurements of the rates for B — » plvg were published subs equently tl^l 
This year, results based on the entire CLEO data (9.7 x 10 6 BB pairs) were reported EDI including the 
measurement of the branching ratio for B + — > rjt + vp (charge conjugation average is implied): 

B(B° -> n-l+vt) = (1.33 ± 0.18 ± 0.11 ± 0.01 ± 0.07) x 10~ 4 , 
B(B° -> p-£ + v e ) = (2.17 ± OMtoil ± 0.41 ± 0.01) x 10" 4 , (83) 
B(B+ -> r)£ + u e ) = (0.84 ± 0.31 ± 0.16 ± 0.09) x 10" 4 , 

where the errors are statistical, experimental systematic, form factor uncertainties in the signal, and form 
factor uncertainties in the cross-feed modes, respectively. Rough measurement of the q 2 distributions 
dT(q 2 )/dq 2 for the {nlvi) and (ptvg) modes by splitting the data in three g 2 -bins were also reported. 
BABAR has also measured the decay B — * plvz I I 

B(B° -> p-l+vt) = (3.29 ± 0.42 ± 0.47 ± 0.60) x 10" 4 , (84) 

where the errors are statistical, systematic and theoretical, respectively. 

Extracting \V u b\ from these measurements is done by using quark models, QCD sum rules, and 
quenched Lattice-QCD calculations for the form factors. We discuss the two main contenders, Lattice 
QCD and QCD sum rules, for the semileptonic decays B — > nive, as this involves (neglecting the lepton 
mass) only one form factor, F + (q 2 ), defined as follows: 

(7r(p w )|67 p g|B( P fl)) - ((pb - ^'^ q^ F + {q 2 ) + w| ~ ^ F (q 2 ) (85) 



with F+(0) = F Q (0). 
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QCD sum rules E2I j n particular Light cone QCD sum rules (LCSRs) | 143 | 144 | have been used 
extensively to study the form factors in B — > it transitions (and other related processes) In the 
LCSR, one calculates a correlation function involving the weak current and an interpolating current with 
the quantum numbers of the B meson, sandwiched between the vacuum ((0|) and a pion state (|7r)): 

F„(p, q)=i[ dxe lqx {-K{p)\T{u{x) llL b{x), m b b(0)i l5 d(0)}\0) . (86) 

For large negative virtualities (g 2 ,p.g) of these currents, the correlation function (CF) in the coordinate 
space is dominated by the dynamics at distances near the light cone, allowing a light-cone expansion of 
the CF (hence, the name). In essence, LCSRs are based on the factorization property of the CF into 
non-perturbative light-cone distribution amplitudes (LCD As) of the pion, called <j)™(u,fj,), where u is 
the fractional momentum of the quark in the pion and /j. is a factorization scale, and process-dependent 
hard (perturbative QCD) amplitudes Tjj(fj,, nib, u, q 2 ), where q 2 is the virtuality of the weak current. 
Schematically, the coefficients in front of the Lorentz structures in the decomposition of the CF Ij86|l can 
be written as: 

C(ji, q 2 ,m b ) ~ T h(^ m b, u, q 2 ) ® <f> n (u, fj) , (87) 

n 

where the sum runs over contributions with increasing twist, with twist-2 being the lowest, and the symbol 
<g) implies an integration over the variable u. The amplitudes Tjj have an expansion in perturbative QCD 
(i.e., a s ). The same correlation function can also be written as a dispersion relation in the virtuality of 
the current coupled to the B-meson. Equating the two, using quark-hadron duality, and separating the 
i?-meson contribution from higher excited and continuum states, results in the LCSR. As an illustration, 
the LCSR for the form factor F + {q 2 ) in the lowest order in a s and leading-twist has the formal 

1 _ m l f 1 du m ^-p 2 i 1 -") 

F +(q 2 ) = o 2 t e ~^™ 2 b U / — e Mu, m,) , (88) 

2m B fB J a u 

where / w = 132 MeV, M is a Borel parameter, characterizing the off-shellness of the 6-quark, and (f> n (u, fib) 
is the twist-2 LCDA of the pion 

( t>«(u 1 fi) = Qu{l-u)(l + al{fi)Cl ,2 {2u-l) + ...) . (89) 

Here, C* 2 3/2 (x) is a Gegenbauer polynomial and (/i) is a non-perturbative coefficient (the second Gegen- 
bauer moment) to be determined, for example, from the data on the electromagnetic form factor of the 
pion. The lower integration limit denoted by A is defined through A = (m 2 ~ q 2 )/{s^ — q 2 ), where Sq 
is determined by the subtraction point of the excited resonances and continuum states contributing to 
the dispersion integral in the B channel. Assuming quark-hadron duality, this subtraction is performed 
at (p + q) 2 > sf. 

In principle, given the assumption of quark-hadron duality, the LCSRs can be made arbitrarily 
accurate, by calculating enough perturbative and non-leading twist contributions to the CF. In practice, 
this framework has a number of parameters (such as fi, M, df(fi), sjf , mj), whose i mprecise kno wledge 
restricts the precision on the CF. Typically, the state-of-the-art LCSRs for F + (q 2 ) D-45 | 146 | 1 ^3 which 
include a s corrections to the leading- twist and part of the twist-three contributions, and tree level for 
rest of the twist-three and twist-four, have an uncertainty of ~ =b209o IS ^ | 148 | anc j y. j g p ro bably difficult 
to make these estimates more precise. 

In Lattice QCD, one calculates a three-point correlation function involving interpolating operators 
for the B and tt mesons and the vector current ~ ^bl^q- In the limit of a large time separation, the 
correlation function has the following behaviour 

C„(p, *, t f , t., U) = Zl' 2 {B J^^ } e~ E ^ e-z^ , (90) 
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where Eb{E^) is the energy of a B(w) meson with the three-momentum k(p) and Z'g 2 (Z^ 2 ) is the 
external line factor calculated from the two-point correlation functions involving the interpolating B{tt) 
fields. Since, one has to go to large time separations to suppress the continuum contribution, one is forced 
to restrict the pion momentum in B — > itlvi decay to low values. Hence, in lattice calculations, there 
is an upper limit on this momentum, |p max |, prescribed by the requirement to keep the statistical and 
discretization errors small. There is also a lower limit on \p\ dictated by the difficulty in extrapolations 
in \p\ and li ght q uark masses. Thus, for example, the FNAL Lattice QCD calculations for the B — > 7r 
form factors have as cut-offs |p max | = 1.0 GeV and |p min | = 0.4 GeV. This translates into a limited 
q 2 range q^ ain < q 2 < g^ ax close to the zero-recoil point. 

Recalling that the differential decay rate for B — > irlve. is given by 

g (J? ^ )s «^, (91) 
dp 2An A E 

where E = p n .pB /fij is the pion energ y in t he B- meson rest frame, one calculates the dynamical part 
on the lattice over a limited region of \p\ ESI Defining 

/•|p max l T) i \F+(E')\ 2 
T B (\p™%\p™\) = J ^ dp P 1 g )l , (92) 

one combines the theoretical rate with the experimental measurements in the same momentum range of 
the pion to arrive at the following relation for | V^,b | 2 , 



, 12tt 3 1 f lP 1 dT(B -> irivi) 

1 T B (\p^-\,\p^\) J„ dp Tp ' 



This avoids the need to extrapolate to higher pion momenta (or low q 2 ) . The practical problem in using 
(I93[l is the paucity of experimental data in low |p|-region, as the differential decay rate has a kinematic 
suppression for low pion momenta. 

Alternatively, one has to use models fo r q 2 e xtrapolation of the Lattice results to lower values of q 2 . 
This is done, for example, by the UKQCD HSDl an d the APE collaborations ESI which make use of the 
LCSRs (discussed above) to constrain the form factors at lower values of q 2 . 

Apart from this, there are also other systematic differences among the various Lattice calculations 
of the B — > 7r form factors, the most important of which is related to the fact that for the current lattice 
spacing a one has m^a > 1. To control lattice spacing effects, one has to do the calculations for values 
of the heavy quark mass toq much smaller than nib and then extrapolate from tuq ~ 1 — 2 GeV, where 
the lattice data is available, to the 6-quark mass. This, for example, is done by th e UKQCD l ^0 | and 
APE ESI collaborations. A different route is taken by the FNAL Lattice group 0321 i n which HQET is 
applied directly to the Lattice observables, using the same Wilson action for fermions as adopted by the 
other groups, but adjusting the couplings in the action and the normalization of the currents, so that 
the leading and the next-to- leading terms in HQET are correct. This allows to perform the calculations 
directly at ttiq — mj. However, there are still some open issues in this approach what concerns the 
non-perturbative matching of the lattice with the continuum. Finally, JLQCD USB a l so uses the Fermilab 
NRQCD approach. The result of the four Lattice groups for the B — > ix form factors, F+(q 2 ) and Fo(q 2 ), 
are shown in Fig. HJ1 taken from the review by Becirevic HS2. These calculations agree at the level of 
±20%, though this consistency is less marked for the form factor Fo(q 2 ). 

The results of the Lattice groups (APeESI UKQCD USB, FNAL Oil and JL QCD ^SH) have also 
been fitted by a phenomenological form for F + ^(q 2 ) due to Becirevic and Kaidalovl ^ 1 

F (n 2, = no) . , _ f(o) 

+ [q ' {l-?/m%.)(l-a ir (P/m%.)' ' I - q 2 / (^m 2 B ,) > ^ 

involving three parameters F(Q), and f3„. The best-fit solution for F + (q 2 ) and F (q 2 ) is shown by the 
dashed curve in Fig. This figure also shows the prediction obtained by the LCSRs 12=1 The resulting 
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Figure 6. Summary of the current (unquenched) Lattice calculations of the B — » n form factors Fo(q 2 ) and F+(q 2 ) and 
from the Light Cone QCD sum rules. The dashed curve shows a parametrization by the Becirevic-Kaidalov model. (Figure 
taken from Becirevicl " I.) 



theoretical description for the form factors from lattice QCD and LCSRs is strikingly consistent, albeit 
not better than ±20%. Further details can be seen elsewhere l ^^ l 

In future, with more CPU power at their disposal, it should be possible to increase the pion momentum 
range accessible on the lattice, allowing a larger and statistically improved overlap of the lattice results 
with the experimental data on B ^ nivg. Also, by using the FNAL method of applying HQET directly 
to the Lattice observables, or else by doing simulations at larger values of uiq than is the case right now, 
theoretical errors on the FFs can be reduced to an acceptable level in the quenched ap prox imation. There 
are also other techniques being developed for treating heavy quarks on the Lattice EMI. The last step 
in this theoretical precision study will come with the estimates of unquenching effects using dynamical 
fermions; first unquenched results for B — ► 7r form factors are expected soon. 

After this longish detour of the currently used theoretic al fra mework, we return to the extraction of 
\V u b\ from exclusive semileptonic decays. The CLEO result EDI quoted below is obtained by using the 
quenched Lattice QCD results for q 2 > 16 GeV 2 and the LC QCD sum rules for lower values of q 2 , 

\V ub \ = (3.17 ± 0.ir±g:if±g*J| ± 0.03) x 10" 3 [CLEO (exclusive)] . (95) 

BABAR uses the theoretical decay width calculated using Lattice QCD, LC QCD sum rules, and 
three quark models to estimate the form factors. The combined result is the weighted average of these the- 
oretical approaches, where the weight is obtained by the theoretical uncertainty, and an overall th eoret ical 
uncertainty is assigned by taking it to be half of the full spread over these models. The result is ^H 

\V ub \ = (3.64 ± 0.22 ± 0.25±g;||) x 10~ 3 [BABAR] . (96) 

The two measurements (|95ll and (19611 arc consistent with each other, and they have been averaged by 
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Schubert ^2 to yield a value of \V u b\ from the exclusive decays, 

iKblcxci = (3.40±°;24 ± 0.40) x 10~ 3 . (97) 

However, this value lies below \V u b\ measured from the inclusive decays B — > X u £i>£, whose current 
average is given in l|82|) . The mismatch in the values of \V u b\ from the inclusive and exclusive decays 
is roughly about 20% and has to be resolved as more precise data and theory become available. The 
possibility that in the quenched Lattice QCD and LCSR estimates, the form factor F + (q 2 ) is estimated 
too high by about 20% can not be excluded at present. 

Digressing from the discussion of \V U b\, we remark that this trend is also seen in the comparison of 
data on B — > K*^ with the LC-QCD sum rule estimates of the B — > K* form factors. To put this in 
a quantitative perspective, we recall that the current branching ratios for the B — > K*^ decays are^^ 
(again charge conjugated averages are implied) 

B{B a -> K* "/) = (4.17 ±0.23) x 10~ 5 , B(B~ -> K*~i) = (4.18 ± 0.32) x 10~ 5 . (98) 

The corresponding theoretic al rat es have been calculated in theNLO accuracy | 155 | 156J 15_7] using the 
QCD- factorization framework ^ ^ 1 An updated analysis based on ^^^ l (neglecting a small isospin violation 
in the decay widths) yields 

where the default value for the form factor T 1 K *(0, nib) is taken from the LC-QCD sum rules EH1 b ; 
and the pole mass po i is the one-loop corrected central value obtained from the MS 6-quark mass 
fhb(mb) = (4.26±0.15±0.15) GeV in the PDG reviews^. Since the inclusive branching ratio for B — > X s ^ 
in the SM agrees well with the current measurements of the same (discussed below), the mismatch in the 
estimates of the exclusive branching ratios in l|99|l and current measurements (|98[) in all likelihood has 
a QCD origin. Of the possible suspects, form factor is probably the most vulnerable link in the chain 
of arguments leading to (|99|l . Interpreting the factorization-based QCD estimates and the data on their 
face value, good agreement between the two requires T^"(0,fh b ) ~ 0.27 ± 0.02. This is shown in Fig.0 
where the ratio 

R Wx«) = l$ZxJy (100) 

is plotted as a function of Tj (0,fhb). The horizontal bands show the current experimental value for 
this quantity R(K*~//X s ~f) = 0.125 ± 0.015. The allowed values of (0, fh b ) are about 25% below the 
current estimates of the same from the LC-QCD approach (= 0.38 ± 0.05). There is a need to do an 
improved calculation of this (and related) form factors. Along this direction ? SU (3)-breaking effects in the 
K and K* LCDA's have been recently re-estimated by Ball and Boglione l -^vU I This modifies the input 
value for the Gegenbauer coefficients in the .KT*-LCDA and the contribution of the so-called hard spectator 
diagrams in the decay amplitude for B — » K*"f is reduced, decreasing in turn the branching ratio by about 
7% I I The effec t of this correction on the form factor (0,m;,), as well as of some other technical 
improvements^^, has n ot yet been worked out. Updated calculations of this form factor on the lattice 
are also under way HS2l with preliminary results yielding values for (0,mf,) ~ 0.27, as suggested 
by the analysis in Fig. and considerably s maller than the ones from the earlier lattice-constrained 
parameterizations by the UKQCD collaboration Theoretical estimates of the form factors are still in 
a state of flux. Phenomenologically, smaller values of the form factors in B — > tt and B — > p transitions are 
preferred by the data, bringing \ V u b\ from the exclusive decays more in line with the value of this matrix 



''The decay rates in this approach depend on the effective theory parameter, called £ K (0), which is re lated by an 0(a a ) 
relation to the B -> K* form factor by Tf*(0,fh b ) ~ C K * (fh b )£ K * (0), with C K '{rh b ) = 1.05 - 1.08 GS2| To keep the 
discussion simple, we used this relation to express the rates in Tj ff *(0, fh b ). 
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Figure 7. The ratio of the branching ratios defined in Eq. I1UUI . plotted as a function of the QCD form factor (0, mj,) 
and the current experimental measurement of this ratio. The solid lines are the centr al ex perimental and theoretically 
predicted values and the dotted lines delimit the ±1<t bands. (Figure updated from Ref. USSj) 



element measured from the inclusive decays. Smaller value of the B — > K* form factor would also improve 
the agreement between the QCD-factorization based estimates for B(B — ► K*j) and experiments. 

A robust average of | V ub | based on current measurements is expressly needed to deter min e one of the 
sides (R b ) of the unitarity triangle precisely. This is, however, not yet provided by HFAG^l A bonafide 
average is difficult to undertake, as the common (and experiment-specific) correlate d syst ematic errors 
are not at hand, as stated in some of the recent experimental reviews on this subject P2 | 53 | j n an y caS g ; 
the dominant errors on \V u b\ are theoretical. Typically, theory-related error from the inclusive measure- 
ments is of 0(15%) at present (and somewhat higher from the exclusive decays), in comparison with the 
experimental error (statistics and detector systematics) on this q uant ity, which is typically of 0(7%). 
This is reflected in the world averages for \V ub \ presented by Stone and Schubert S2, respectively, 

\V ub \ = (3.90 ± 0.16(exp) ± 0.53(theo)) x 1(T 3 , 

(101) 

\V ub \ = (3.801^ (exp) ± 0.45(theo)) x 10~ 3 . 

Adding the errors in quadrature, the first of these leads to \V u b\ = (3.90 ± 0.55) x 10~ 3 , yielding 
^|Kib|/|Kib| = 14%, and a very similar range if one uses the value given in the second. Thus, the 
matrix element \V ub \ is still considerably uncertain, and we trust that the B factory experiments and 
theoretical developments will make a major contribution here, pushing the error down to its theoretical 
limit 0(5%), mentioned earlier. 

Using the current averages, |V c6 | = (41.2 ± 2.0) x 10~ 3 and \V ub \ = (3.90 ± 0.55) x 10~ 3 , we get 

= 0.095 ± 0.014 =*> R b = \ = 0.42 ± 0.06, (102) 
\V cb \ A \V cb \ 

which determines one side of the unitarity triangle. 
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5 Status of the Third Row of Vckm 



Knowledge about the third row of the CKM matrix Vckm is crucial in quantifying the FCNC transitions 
b — > s and b — > d (as well as s — > d) and to search for physics beyond the SM. The FCNC transitions 
in the SM are generally dominated by the top quark contributions giving rise to the dependence on the 
matrix elements |V^Vt s | (for b — » s transitions) and |V^V{<i| (f° r b — > d transitions). Of these, only the 
matrix element \Vtb\ has been measured by a tree amplitude t — ► Wb at the Tevatron through the ratio 



_ B(t^Wb) = \V tb \ 2 

* b " B(t -> w^) |v id | 2 + |v fs | 2 + |v tb | 2 • 1 UOj 



The current measurements yield El R tb = 0.94+Q;|g, which in turn gives 



\V tb \ = 0.96lg;23 => \Vtb\ > 0.74 (at 95% CL) . (104) 

Thus, this matrix element is consistent with unity, expected from the unitarity relation |V^f,| 2 + 
|V^b| 2 + \Vtb\ 2 = 1, though the current precision on the direct measurement of |Vtb| is rather modest. 
(Unitarity gives \Vtb\ — 0.9992. ) T he precision on |Vt&| will be greatly improved, in particular, at a 
Linear Collider, such as TESLAl ™ \ but the corresponding measurements of |Vf S | and |Vfd| from the tree 
processes are not on the cards. They will have to be determined by (loop) induced processes which we 
discuss below. 



5. 1 Status of | Vtd | 

The current best measurement of \Vtd\ comes from AM Bd , the mass difference between the two mass 
eigenstates of the BS - B®, complex. This has been measured in a number of experiments and is known 
to an accuracy of ~ 1%; the current world average is 1^ AM Bd = 0.502 ± 0.006 (ps) -1 . 

In the SM, AM Bd and its counterpart AM Bb , the mass difference in the B® - B® system, are calculated 
by box diagrams, dominated by the Wt loop. Since (M\y,rn t ) ^> nib, AM Bd is governed by the short- 
distance physics. The expression for AAI Bd taking into account the perturbative-QCD corrections reads 
as follows USa 

AM b, = ^ VB \V td V t * b \ 2 m Bd (f 2 Bd B Bd ) M 2 W So(xt) ■ (105) 

The quantity fj B is the NLL perturbative QCD renormalizati on of the matrix element of the (|AB| = 
2, AQ = 0) four-quark op erato r, whose value is i) B = 0.55±0.01 x t = mf/M^ and So(x t ) = x t ji(xt) 
is an Inami-Lim function EHI with 

, . , 1 9 1 3 1 3 x 2 lna; 

MX) = 4 + 4 (1 = *) 2 (I" ^ 2 (I" ^ • ( 106 ) 

The quantity fi B Bd enters through the hadronic matrix element of the four-quark box operator, de- 
fined as: 

<5°l(Ml - 7 5 )<Z) 2 |B°) = \f%B Bq Ml q , (107) 

with B q = Bd or B s . With AM Bd and f) B known to a very high accuracy, and the current value of the 
top quark mass, defined in the MS scheme, m t (m t ) = (167 ± 5) GeV, known to an accuracy of ~ 3%, 
leading to 8So(xt)/So(xt) — 4.5%, the combined uncertainty on \Vtd\ from all these factors is about 
3%. This is completely negligible in comparison with the current theoretical uncertainty on the matrix 

element f Bd \J B Bd . For example, 0(a s )-improved calculations in the QCD sum rule approach yield E31 

f Bd = (210 ± 19) MeV and EDI f Bd = (206 ± 20) MeV, whereas B Bd in the MS scheme in this approach 
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is estimated as ED Bs d — 1 to within 10%, yielding for the renormalization group invariant quantity 



Bs d — 1.46, and an accuracy of about ±15% on fB d y Bs d 



Lattice calculations of fB d \jBB d are uncertain due to the chiral extrapolation. This is shown in 
Fig.|H|from the JLQCD collaborationEl 

in which the quantity $ f B = f Bq y/Ms q is plotted for q = d,s 
as a function of the pion mass squared, with both axes normalized with the Sommer scale tq determined 
from the heavy quark potential at each sea quark mass. The lattice calculations in this figure are done 
with two flavours of dynamical quarks u and d, for the u and d quark masses in the range (0.7 - 2.9)m s , 
with m s being the strange quark mass. The solid line represents a linear plus quadratic fit in (rom.^) 2 , 
which describes the lattice data well. This fit, however, does not contain the chiral logarithmic term, 
predicted by the Chiral perturbation theory 1 1 



!B d ^fm d 3(l±3 ff 2 ) 



in 



2 



(f Bd ^)^ 1 ~ ln # + "' ' (108) 

where terms regular in m 2 are omitted, f n = 130 MeV, and g is the B* Bit coupling in chiral perturbation 
theory. A recent lattice calculation HZa g ives Qb* b* — 0.58 ± 0.06 ± 0.10. A related quantity go* dtt has 
been determined from D* — > Dtt decay E3I g D * D7r = 0.59 ± 0.01 ± 0.07. The value of g is fixed at 
g = 0.6 in drawing the three curves with the chiral behaviour H108|) with the three values of the hard 
chiral cutoff: /j, = 300 MeV (dotted curve), fj, = 500 MeV (thin dashed curve) and /i = oo (thick dashed 
curve). Lattice data with the currently used high values of the dynamical quarks is not able to verify the 
chiral logarithm. The data are not inconsistent with such a behaviour either. 
The chiral behaviour of the bag constant is given by the expression 

% = 1 _ ('-V) nt < + 
B ( „> 2 (471-/* ) 2 M 

the coefficient (1 — 3g 2 )/2 is numerically small (= —0.04), as opposed to the coefficient in $/ Bd for which 
3(1 + 3.g 2 )/4 = —1.56, with g = 0.6. Hence, extrapolation of the lattice data to the small quark masses 
poses no problems for the bag parameters. 

Takin g thi s into account, the unquenched lattice QCD calculation from the JLQCD Collabora- 
tion yields LLLzl 



f Bd ^B Bd - 215(11)(3 3 )(15) MeV, (110) 

where the first error is statistical, the second (asymmetric) is the uncertainty from the chiral extrapolation 
and the last is the systematic error from finite lattice spacing. The largest error is from the chiral 
extrapolation in fg d , which in the conservative estimate of the JLQCD collaboration could be as lar ge as 
-10%, with f Bd = 191(10)(±g 9 )(12) MeV. For further discussion, see the recent reviews by Becirevic* \ 
Kronfeld GIDand WittigESl 

We shall use the unquenched lattice result in l|110fl by adding the errors in quadrature and sym- 
metrizing the errors, getting 



! Bd \jB Bd = 210 ±24 MeV. (Ill) 

The dependence of \V t d\ on the various input parameters can be expressed through the following numer- 
ical formula 



\V td \ = 8.5 x 10" 
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Br 



AM 
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Figure 8. Chiral extrapolation of ^f Bd (filled circles) and &f B (open squares). The quadratic extrapolation is shown by 
solid lines, while the fits with the hard cutoff chiral logarithm are shown for fi = 300 MeV (dotted curve), 500 MeV (thin 
dashed curve ) and oo (thick dashed curve). Quenched results are also shown (triangles). (Figure taken from the JLQCD 
collaboration LLZal. ) 



where the d efault value of So(xt) corresponds to m t (m t ) = 167 GeV, and the dependence of this function 
on m t is So(xt) ~ x]- 52 . To get the ±ler range for \V t d\, we vary the input parameters within their 
respective ±lcr range and add the errors in quadrature. This exercise yields 

\V td \ = (8.5 ± 1.0) x 10~ 3 . (113) 

In future, unquenched lattice results for fg d and other quantities will be available at smaller values of 
the dynamical quark masses (than is the case in the current JLQCD calculation), allowing to check the 
chiral logarithmic behaviour of fs d , or at least reduce the error associated with this extrapolation. 

Knowing \ Vtd\ from 1113(1 . and \V c b\, and A from previous sections, one can determine the other side 
of the UT, which has the following central value: 



A \ Vcb I 



|Vm| 



8.5 x 10 



-3 



0.041 



V, 



(114) 



Taking into account the errors (and taking symmetric errors on |Vtd|), we get Rt = 0.93 ± 0.12. 
5-2 \Vtd\ from B — ► pj decays 

Independent information on \Vtd\ (more precisely on p and ff) will soon be available from the radiative 
decays B — * (p, u>)~/. There is quite a lot of theoretical interest lately in this process, starting from 
the earlier papers a decade ago ^^ ^^, where the potential impact of these decays on the CKM phe- 
nomenology was first worked out using the leading order estimates for t he pen guin amplitudes. Since 
then, annihilation contributions have been estimated in a number of paper ™^^ ^ , and the next-to- 
leading order corrections to the decay amplitudes have also been calculated ESESI Deviations from the 
SM estimates in the branching ratios, isospin-violating asymmetry A ±0 and CP-violating asymmetries 
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-4cp(p ± 7) an d -4cp(/0°7) have also been worked out in a number of theoretical scenarios | 18 4 | 185 | 186 | 
These CKM-suppressed radiative penguin decays were searched for b y the CLEO colla borat ion U^ffi and 
the searches have been set forth at the B factory experiments BELLE^EI anc j BABAR ^^ I The current 
upper limits at 90% C.L. (averaged over the charge conjugated modes) are given in TableQ] 



Table 1. 90% c onfide nce leve l upp er limits on the branching ratios (in units of 10 6 ) for the decays B 
from the CLEoEHB BELLE Esfil an d BABArCHI collaborations. 



fry and B 





B(B+ -> p+ 7 ) 


B{B° -> p° 7 ) 


B(B Q -> ury) 


CLEO (9.1 fb- 1 ) 


13.0 


17.0 


9.2 


BELLE (78 ftr 1 ) 


2.7 


2.6 


4.4 


BABAR (78 fb" 1 ) 


2.1 


1.2 


1.0 



The BABAR upper limits on £(5+ -> p + 7 ) and B(B° - 
symmetry to yield an improved upper limit I I 

B(B -> p 7 ) < 1.9 x 10" 



p 7) have been combined using isospin 



(115) 



Together with the current measurements of the branching ratios for B — > K**y decays, studie d ear lier, 
this yields a 90% C.L. upper limit on the isospin-weighted and charge-conjugate averaged ratio USall 



B(B -» P7 ) 



< 0.047 



(116) 



B{B -> K*-y) 

The branching ratios for B — » P7 ha ve been calculated in the SM at next-to- leading order E 5 | 156 | 
in the QCD factorization framework ^ SBI . As the absolute values of the form factors in this decay and in 
B — > i^* 7 decays discussed earlier are quite uncertain, it is advisable to calculate instead the ratios of 
the branching ratios 



i? ± (p 7 /^*7) 



B(B d 



B(B± - 
B{B° 



K*± 7 )- 



B(B" -> *f*° 7 )' 



(117) 



(118) 



The results in the NLO accuracy can be expressed as 



i? ± (p 7 /^*7) 



Id 



(M 2 B 



2\3 



M 2 ) 



(M 



R°(m/K*i) = 2 



2 (Mi 



(Mi 



C 2 (l + AR (e° A ,p,fj)) 



(119) 



where C = Tf (0) /T^* (0) , with Tf (0) and * (0) being the form factors evaluated at q 2 = in the decays 
B — ► pT and £? — * ^* 7 , respectively. The functions Ai? ± (e^ : , p, 77) and Ai?°(e^, p, 77), appearing on the 
r.h.s. of the above equations encode both the 0(a s ) and ann ihilation contributions, and they have a non- 
trivial dependence on the CKM parameters p and Updating them, incorporating also a shift 
in the quantity called A^ 1 , related to an integral over the i?-meson LCDA, A^ 1 = J^dk/k cp + (k,n), 
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which has been evaluated in the QCD sum rule approach rece ntly b y Braun and Korchemsky EDI 
= (2.15 ± 0.50) GeV-\ the result for the functions in (EH isG&H 



A 



AR ± = 0.056 ± 0.10 , AR Q = -0.010 ± 0.064 , 



(120) 



where the uncertainties reflect also the variations in the CKM parameters p and rj, for which the ranges 
p = 0.21 ± 0.09 and fj = 0.34 ± 0.05 have been used. 

Theoretical uncertainty in the evaluation of the ratios R ± (pj/K*^f) and R°(p^/K*j) is dominated 
by the imprecise knowledge of the quantity £. In the SU(3) limit £ = 1; SU(3)-breaking corrections have 
been calculated in several approaches, including the QCD sum rules and Lattice QCD. In the earlie r 
calculations of the ratios, the follo wing ranges were used C = 0.76 ± 0.06 (by Ali and Parkhomenko^^J), 
C = 0.76 ± 0.10 (Ali and LunghiE^ and 1/C = 1.33 ± 0.13, leading to C = 0.75 ± 0.07 (Bosch and 
BuchallaD3S3). These ranges reflect the earlier estimates of this quantity in the QCD sum rule approach 
| 182 | 192 | 193 | 194 | an( j i nc Ji ca tes substantial SU(3) breaking in th e B — > V form factors. Now there exists an 
improved Lattice estimate of this quantity, with the result^^^ = 0.9±0.1, which is within la compatible 
with no SU(3)-breaking! We conclude that £ is at present poorly determined. It is essential to calculate 
it precisely if the measurement of R(pj/ K*j) is to make an impact on the CKM phenomenology. 

To illustrate the impact of the current bound on R(pj/ K*j), we use the following estimate 



C = 0.85 ±0.10. 

Including the uncertainties from other input parameters, the updated results are ^^ 

R ± {p 1 /K* 1 ) = 0.033 ± 0.012 , 

R°(pj/K*j) = 0.016 ± 0.006 . 

Combining these with the measured values of the branching ratios B{B ± — > K* ±n /) and B(B° 
the predictions for B(B ± — > p ± -f) and B(B° — > p°-f) are as follows: 

p± 7 ) = (1.36 ± 0.49[th] ± 0.10[exp]) x 10~ 6 , 



(121) 
(122) 



(123) 



B(B ± 

B(B° -> p°7) = (0.64 ± 0.23 [th] ± 0.04[exp]) x 10~ 6 , 

and B(B° — ► p°j) = B(B° — > CJ7) for the stated range of theoretical uncertainty. Comparing these 
predictions with the present experimental bounds given in Table ^ we expect that all these branching 
ratios lie within a factor 2 of the current experimental bounds, and hence will be measured soon. 

The isospin-weighted and charge-conjugation averaged ratio R^j/K*^) is given by the following 
expression in the SM^I 

A 2 C 2 



R( P1 /K* 7 ) 



^ -> {Ml d -mlf 



[1 



(M^-m^) 
Re lG x {p,e±)+r, 2 G 2 (e± 



4p] 2 + (i-4)V 



+ (4 



4) 



(124) 



where the NLO contribution are introduced through the functions: 
2 



Gi(p,e A ) 



(O)cff 



{(1 - pf [a%>> - Ag*'] + p(l - p) [A- + EA AW*] + p 2 e A A u ) 



G,(e A ) = -^{A^-A^' -A* + 



A u _ A (l)t 



} 



(125) 



Here, ca represents the annihilation contribution, estimated as ~ 0.3 ± 0.07lS2 w ith < due to 
being colour- and electric charge suppressed, and the other quantities in l|124|) and (|125|) can be seen in 
the literature E2I 
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The dependence of the ratio R(frf/K*j) on |Vtd|/|V ts | is shown in Fig. |5J Note, that this deviates 
from a quadratic dependence, which holds only if one neglects the annihilation and 0{a s ) corrections. 
Including these corrections, we have given the dependence of R(pj/ K*j) on p and fj explicitly. The solid 
curve corresponds to the central values of the input parameters, and the dashed curves are obtained by 
taking into account the ±1ct errors on the individual input parameters in R(p r f / K*j) and adding the 
errors in quadrature. The current experimental upper limit on R(pj/ K*j), given in l|116|) . is also shown. 
Taking the least restrictive of the three theoretical curves, the current experimental upper limit yields 
iVy/lVy < 0.28, to be compared with the SM range |V t(i |/|V" ts | = 0.21 ± 0.03, shown as the region 
between the two dotted horizontal lines. Conversely, constraining the ratio | V^^ | / 1 | in the SM range, 
we get R{frf/K*i) = 0.032 ± 0.008. 
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Figure 9. Ratio of the branching ratios for the decays B — » p7 and B — > K*^, R(py / K* -y) , plotted as a function of \Vtd/Vt a \. 
The current experimental upper limit R(p~f / K*~f) < 0.47 is shown as the horizontal line. The vertical dotted lines demarcate 
the ±lc range of |Vtd/Vt s | from the SM-based unitarity fits. The solid curve corresponds to the central values of the input 
parameters and the dotte d cur ves delimit the ±lcr range in the NLO corrected calculations for R{py/K*y) in the SM. 
(Figure updated from Ref.li^l) 



5.3 Present status of \Vt s \ 

There are two measurements at present which yield indirect information on | Vt s | , namely the lower bound 
on the B® - B® mass difference AMb s and the measurement of the branching ratio B(B — ► X s j). We 
discuss them in turn below. 

The expression for AMb b in the SM can be obtained from the one for AMs d in l|105|) by the 
replacements: E>d — > B s , V t d — >• V ts - However, as opposed to AM^, AMb s is not yet measured, and the 
current lower bound (at 95% C.L.) is ^ AMb s > 14.4 (ps) _1 . This can be converted into a lower bound 

on | Vt s | , knowing fg t \J ' Bb s - The unquenched lattice QCD calculation from the JLQCD Collaboration 
for f Bs and f Bs \[~B~b~ s yields E2 

= 215(9)(±g)(13)(±g) MeV, 

" (126) 
fB.\f&B.= 255(10) (11) (17) C±S) MeV, 
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where the errors have the same origin as in (|110|) , and the additional error here is due to the ambiguity in 
the determination of the s-quark mass. An average of the JLQCD, MILC and CP-PACS data gives' 



f Bs \J Bb s = 254 ± 13 ± 14 ± 13 MeV. However, a recent calculation of the coupling constants J Bb and fo s 
in the unqucnchcd lattice QCD including the effects of one strange sea quark and two light sea quarks 
by Wingate et aZ. ESI yields, 

f B = 260 ± 7 ± 26 ± 8 ± 5 MeV , 

(127) 

f Ds = 290 ± 20 ± 29 ± 29 ± 6 MeV , 

where the errors are respectively due to statistics and fitting, perturbation theory, relativistic corrections, 
and discretization effects. The result for f Bs in (|127fl is typically 20% higher compared to the one from the 
JLQCD collaboration given in Ijl26|l . As both of these calculations are based on the NRQCD framework 
for heavy quarks, the difference between the two lies in the details of the lattice simulations, such as 
the dynamical quark masses used and n/, which are different in the two approaches. Based on these 



comparisons, we conclude that the current lattice precision on f B , y B B , is of order 20%. 

It has become cus toma ry to use the ratio of the mass differences AM Bd / AM Bs to constrain |Vtd|/|Vfs| 
from the SM relational 

AM Bs _ M Bs \V*V ts \ 2 

AM Bd ? M Bd \V t lV td ? ' 1 1 

where 



£ = (129) 

fB d V B Bd 



Theoretical uncertainty in £ is arguably smaller compared to the one in f Bs y Bb s , as in the SU(3) limit 
£ = 1, and the uncertainty is really in the estimate of SU(3)-breaking corrections. Thus, <5£ ~ 10% is not 
an unrealistic error on £. Current estimates can be exemplified by the unquenched lattice calculations of 
£ from the JLQCD collaboration 

£ = 1.14(3)(±; 3 )(2)(tg). (130) 
The single largest uncertainty in £ is due to the chiral extrapolation - the same source as in the estimates 



°f fB d \j Bb a - Symmetrizing the errors, JLQCD result yields 

£ = 1.19 ±0.09. (131) 

It should b e rem arked that the quantity S = £.§ ■*/ Ik is useful to estimate £, as the chiral logs largely 
cancel in 5 Using the JLQCD data tL7JI, Kronfeld quotes £ = 1.23 ± 0.06, with the error 

increasing if one includes the preliminary HPQCD results^, yielding £ = 1.25±0.10. The errors in this 
and the one in H131|l are almost the same and about 8%. Thus, £ has a value in the range 1.1 < £ < 1.3. 

The constraint that a measurement of (or equivalently a bound on) AMs a provides on R t can be 
expressed as follows 



Rt = 0.90 



£ 
1.20 



^J^i, (132) 
AM Bs A/ 0.50/ps ' v ; 



where the default value of AAI Bs is the best-fit value from the CKM unitarity fits, discussed in detail 
later c . The present bound AM Bs > 14.4 (ps) -1 yields 

\Vtd\ 



\V ti 



>0.22 => \V* b V ts \ > 0.034, (133) 



C I acknowledge the help provided by Enrico Lunghi and Alexander Parkhomenko in updating the CKM unitarity fits. 
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where the last inequality follows from using |V^|Vtd| > 7.5 x 10 3 , determined earlier. This is to be 
compared with the CKM-unitarity constraint | Vtf> V^* | = (V^V^I + 0(X 2 ), which for the central value of 
\Vcb\ = 0.041, predicts \V ts \ = 0.04. 

5.4 Determination of \Vt s \ from T(B — > X s j) 

We now discuss the determination of |Vt s | from B(B — ► X s j). The effective Hamiltonian which governs 
the transition B — > X s j in the SM is 



H eS = -^V: s V tb Y,C i {n)O i 

v l—l 



(A*), (134) 



which is obtained by integrating out all the part icles that are much heavier than the 6-quark. The 
operators Oi can be seen, for example, in Ref. LLSua and unitarity of the CKM matrix is used to factorize 
the CKM-dependence of H e g in the multiplicative product V t * s V t b- 

The decay rate for B — > X s j is calculated by taking into account the QCD perturbative and power 
corrections. What concerns the perturbative corrections, there are two effects: (i) renormalization of the 
Wilson coefficients d(Mw) — ► £*(/■*&)> where fib ~ 0(m&), and (ii) perturbative corrections to th e matrix 
elements of the operators (Oi). From step (i), one has a perturbative expansion for C,; (^&) EH01 

CM = cfV) + >0*) + • • • • (135) 

In the leading order, i.e., without the QCD corrections, (/itf,) = Ci(Mw) and C- (/x;,) = 0, and of 
the Wilson coefficients Ci(Mw) only Cr(Mw), corresponding to the electromagnetic penguin operator 
07 = errib/ gs(sLO~ til 'bii)F fl „, is relevant for the decay b — ► sj. As Cj(Mw) is dominated by the (virtual) 
top quark contribution, the amplitude M(b — * sj) is proportional to X t = VtbV t *. In this order, the 
contributions from the intermediate up and charm quarks are negligible, being power suppressed due to 
the GIM mechanism I3J So, in the leading order, the decay width T(B — > X s j) depends quadratically 
on \V*V tb \. 

However, including QCD corrections, the power-like GIM-suppression of the intermediate up and 
charm quarks is no longer operative. The reason for this is that QCD forces a very significant operator- 
mixing between the operator O2 = (slT/jCl^Ci^^&l), whose Wilson coefficient C'2 is of order 1, and the 
electromagnetic penguin operator O7, whose Wilson coefficient is much smaller, Ct(ijl) <C 1. Hence the 
contribution from the intermediate charm state becomes numerically very important bringing with it the 
dependence of the decay amplitude on A c = V c bV* s . The contribution from the intermediate u-quark 
can always be expressed in terms of Xt and A c , using the unitarity relation T, UtCtt Xi = 0. However, on 
noting that X u /Xt — 0(A 2 ), A M can be dropped, to an excellent approximation, yielding A c = — At, and 
the electromagnetic penguin amplitude b — > 57 factorizes in At (= — A c ). 

Following this line of argument, one fixes the value of Xt = — X c = (41.0 ± 2.0) x 10~ 3 in calculating 
the SM decay rate for B — > Xsj. Thus, for example, in the MS scheme, taking into account the next- 
to-lead ing orde r (in a s ) and leading order (in 1/m 2 . and 1/m 2 ) power corrections, the SM branching 
ratio is 76 201 : 

B{B^ X sl ) = (3.70 ±0.30) x 10" 4 , (136) 
to be compared with the current experimental world average (based on the CLEO, ALEPH, BELLE and 



BABAR measurements) of this quantity 

B(B -> X sl ) = (3.34 ±0.38) x 10~ 4 . (137) 
The consistency of the two implies that the CKM unitarity, implemented through the unitarity relation 

V ub V* s + V cb V c * s + V tb V t * s = , (138) 
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holds within experimental and theoretical precision. Note, that this is a different unitarity relation than 
the one given in (|47|1 and shown in Fig. [5J While it does not provide any information on the CKM 
parameters p and fj, or for that matter on a, /3 and 7, it involves the CKM matrix element Vt s , apart 
from the other known ones, and its best use is to determine \Vt s \ and its argument 8~/ s . 

We address the question how to convert the information on B(B — ► X s ^f) to determine |Vt a |. To do 
that, we have to keep the individual contributions from the intermediate u, c and t quarks with their 
respective CKM dependencies, X u , X c and At, and not invoke unitarity, in calculating the decay width 
T(B — > Xsj). This is just a different book keeping of the partial contributions to the decay amplitude in 
the SM. Dropping small numerical contributions (< 2.5%), current measurements of B(B — > X s ^f) yield 
the following relation 

|1.69A U + 1.60A C + 0.60A t | = (0.94 ± 0.07)|F c6 | . (139) 

Note the much larger coefficient of A c compared to the coefficient of At, reflecting the large O2 - O7 mixing 
under QCD renormalization. Note also that the relative signs of A c and A t are opposite, which means 
that a destructive interference between the charm and top quark contributions is absolutely essential to 
explain the observed branching ratio for B — > X s j. Solving this equation with the known values of the 
CKM factor s, A e = (41.2 ± 2.0) x 10~ 3 and X u (which is complex) from the discussions earlier (or from 
PDG) yields EOS 

At = V tb V t * = -(47 ± 8) x 10- 3 . (140) 

The reason for the large error on At is, apart from the experimental precision on B(B — > X s "f), the 
relatively small coefficient of this term in l|139|) . Within measurement errors, this determination of X t 
is consistent with the CKM unitarity expectations \Vt s \ = \V c b\ + 0(X 2 ), though it is less precise at 
present than \V c b\. However, due to the interference of the terms proportional to A c and At, B{B — > 
Xg'y) determines the relative sign of At and A c . The sign in (|140|l is in accord with the Wolfenstein 
parametrization given in which has At = — AA 2 . Hence, A is positive definite. 

6 Summary of the Current Status of Vckm an d the Weak Phases 

The current knowledge of the magnitudes of the CKM matrix elements that we have discussed in the 
previous sections is summarized in Table El Note, that these are direct measurements in the sense that 
unitarity has not been used in arriving at these entries. The corresponding values obtained on using the 
unitarity constraints are lot tighter, as also discussed above for some specific matrix elements. 

This information can also be expressed in terms of the Wolfenstein parameters and the sides of the 
unitarity triangle in the SM: 

A = 0.83 ±0.04, 

A = 0.2224 ±0.0020, 

(141) 

Rt = 0.93 ±0.12, 
R b = 0.42 ± 0.06 , 

where the range of Rt given above is coming from AMs d = 0.503 ± 0.006 (ps) _1 . The current lower 
bound on AAfg s also gives a constraint on R t , which already is quite effective in restricting the allowed 
p - fj space in the SM. Finally, the quantity called R(fr//K* r y) also constrains p and fj, but the current 
upper limit is less effective than either AMs d or AMg t . 

So far, we have seen that the unitarity of the CKM matrix, as determined through the magnitudes 
of the matrix elements in each row or each column, holds with deviations which are statistically not 
significant. In the rest of this section, we use the information on the CP-violating asymmetries in the 
Kaon and £?-meson systems to see first the consistency of the data in terms of the sides and the angles of 
the unitarity triangle, and then, as the final step, we undertake a fit of all the relevant data to determine 
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Table 2. Summary of current measurements of \Vij\. 



\Vij\ Value <5|^'l/l%l 



\Vud\ 0.9739 ±0.0005 5 x 10~ 4 

\V US \ 0.2224±0.0020 9 x 10~ 3 

\V ub \ (3.90 ± 0.55) x 10~ 3 14% 

\Vcd\ 0.224 ±0.016 7% 

\V CS \ 0.97 ±0.11 11% 

\V cb \ 0.041 ±0.002 5% 

\V td \ (8.5 ± 1.0) x 10~ 3 12% 

| Vi S | 0.047 ±0.008 17% 

\V tb \ 0.96i°il 20% 



the apex of the unitarity triangle shown in Fig. [5J This will allow us to update the predictions for some 
interesting quantities which have either not been measured yet, such as AMb, , or not precisely enough, 
such as the angles a and 7. 

6. 1 Precise tests of the CKM theory including CP-violating phases 

In addition to the constraints that we have already given in ltlTT|) . (fT3"2"|) (for AM S J, the theoretical 
expression for R(p-f / K*-f) given in (|124|l and l|125|l . and the current experimental bound on this ratio 
in l|116(l . there are three precise measurements involving CP violation in the K- and B- meson sectors 
providing constraints on the CKM parameters. We discuss them briefly here. 

The observed CP-Violation in the Kl — > mr and K$ — + tttt decays yield the following information on 
the quantities \ek\ and Re(e'/e)Gl 

\e K \ = (2.280 ±0.013) x 10~ 3 , 

(142) 

Re(e'/e) = (16.6 ± 1.6) x 10~ 4 . 

The value quoted for Re(e'/e) is the world average from the NA48EEland KTEV EOS collaborations, 
including also the earlier results from their forerunners, NA31 and E731, respectively. While Re(e'/e) is 
a benchmark measurement in flavour physics, as so far this is the only well established example of CP 
violation in decay amplitudes, unfortunately its impact on the CKM phenomenology is muted due to the 
imprecise knowledge of the hadronic quantities needed to extract the information on the CKM parameters 
quantitatively. Given the hadronic uncertainties, which admittedly are not small, the measured value of 
Re(e'/e) is in agreement with the SM estimates. For further details and discussions, the interested reader 
is referred to a recent review on this subject by Buras and JaminEDSl where references to the original 
theoretical papers in the analysis of Re(e'/e) can also be found. 

Recently, time-dependent CP asymmetries in B — ► J/ipKs and related decays have been measured. 
Denoting these asymmetries generically by a^K s {i)i one can measure sin 2(3 (or sin2</>i) from the time 
dependence of the asymmetry 

a^ Ks (t) = a^ Ks s'm(AM Bd t) = sin 2/3 s'm(AM Bd t) ■ (143) 
As opposed to the the oreti cal predictions for \ejc\ and Re(e'/e), the CP asymmetry a^K s {t) is free of 
hadronic uncertainties EH3 which allows to write down the above expression. Current measurements 
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s 



of a^Ks are dominated by the BABArEHI ( ai>Ks = 0.741 ± .0.067 ± 0.033) and BELLeEHI ( a ^ K 
0.733 ± 0.057 ± 0.028) measurements, and the world average E2E3I 

a ipKs = 0.736 i Q.049 , (144) 

is in good agreement with the predictions in the SM^l We shall also quantify this agreement below. 

We now discuss the constraints on the CKM parameters in the SM that follow from the measurements 
of \ex\ and a^Ks- The expression for \ejc\ in the SM, including NLO corrections is l ^^ l 

\e K \ = C K B K (A 2 \ 6 ij) (x c {V3f 3 (x c ,xt) - Vx) + % x t f 2 (x t ) A 2 A 4 (1 - p)j , (145) 

where C K = G 2 F f K m K M^/6V2TT 2 AM K , with El AMk = (3.490 ± 0.006) x 10~ 12 MeV the K° - W 
mass difference and fx — (159.8 ± 1.5) MeV the if -meson coupling constant; Xi — mf/M^, and f 2 (x) 
and fs(x, y) are the Inami-Lim functions EHI f which we have already given f 2 (x) in l|106(l . and fa(x, y) 
is given by the following expression (for x <C y) 

f^'^W^){ 1 + T^^)- (146) 

The rji are perturbative renormalization constants calculated in NLO accuracy l ^^^ l ^l^ l They depend on 
the definition of the quark masses, and the values 772 = 0.57 ± 0.01 and 773 = 0.46 ± 0.05 that will be 
used for numerical analysis below correspond to the definitions wit — TOt(mt) and m c = m c (m c ), which 
represent the quark masses in the MS scheme at their indicated scales. With this definition, the residual 
mi-dependence of r)2, and that of 773 on m t and m c , are negligible, but the dependence of 771 on m c is 
significant. Following Gambino and Misiak in the CERN-CKM proceedings El we shall use 

/ 1 30 \ 1-1 

77! = (1.32 ±0.32) —-— , (147) 
\m c (m c ) J 

and take m c (m c ) = 1.25 ±0.1 G eV. T he quantity Bk is the bag parameter, for which lattice QCD 
estimates yield B K = 0.86(6) (14) EH. 

In working out the constraints from ck numerically, we shall 

take B K = 0.86 ±0.15. 

A useful numerical expres sion showing the constraints that the current value of \ck\ provides on the 
CKM parameters is as follows ESI 

77 [(1 - p) A 2 772 So(x t ) + P c (e)] A 2 B K = 0.187, (148) 

where So(xt) = Xtf 2 (xt) — 2.40 (rat/167 GeV) 152 , and P c {b) summarizes the contribution from the first 
row in 1)145)) . which dep ends on both m c and m t . Taking into account the dependence of r\\ on m c , this 
quantity is estimated asEI2 P c (e) = 0.29 ± 0.07. 

Measurements of a^Ks — sin 2/3 EH3 translate into the following constraints on the Wolfenstein 
parameters 



2(1-/3)71 1 ± Jl - sin 2 2/3 . 

sin 2/3= V or 77= 1 - p ). 149 

(1 — pY + rf sin 2/3 

The constraints resulting from the five quantities [ol^Ksi \ ^k\i AAfe d , AJVfg,, and Rj,) on the CKM pa- 
rameters p and 77 are shown in Fig. 1101 Of these, the allowed bands correspond to ±lcr errors on a^K S i 
|ei£-|, AMB d and and the constraint shown for AMb, is for £ = 1.28, the maximum value in the 
±lcr range given in ()131[1 . with AMb, > 14.4 (ps)^ 1 . Of the two solutions shown for a^K s -, only one 
is compatible with the measurement of Rb, i.e. with the SM, and for this solution we have a consistent 
description of all the data indicated in this figure; the resulting allowed region is shown as a shaded area. 
Note, that this is not a fit, but a simple consistency check of the CKM theory with a large number of 
measurements. The three dotted curves labelled as R(p~f / K*~f) refer to the 90% C.L. bound on this 
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ratio, and show the current theoretical uncertainty in the interpretation of this bound, which is domi- 
nated by the imprecise knowledge of £: C = 0.75 (the leftmost curve), C = 0.85 (the central curve), and 
C = 0.95 (the rightmost curve). 

Fig. 1101 is quite instructive. It shows that the knowledge of Rb is required to distinguish between the 




Figure 10. The constraints resulting from the measurements of a^,n s , the ratio Ri,, AMs d , £k, and the upper bounds on 
AMb, and the isospin-weighted and charged-conjugate averaged ratio R(pj/K*y), in the p - f\ plane. The overlap region 
is indicated by shaded area. 

two allowed solutions for a^K a ■ However, for the solution compatible with the SM, the allowed region in 
the (p, fj) plane is now largely determined by the measurement of a^x s an d AMs d , and the bound on 
AMb s - The constraint from \ck\ is still required and the compatibility of |ejf| and a^K s is an important 
consistency check of the CKM theory. However, as the bound on AMg ( becomes stronger, more so if 
AMb 3 is measured which we anticipate soon, the allowed unitarity triangle could be determined entirely 
from the i?-meson data. This is potentially good news, as the uncertainty on Bk is still quite substantial. 
This figure also reveals that the theoretical uncertainty on the ratio R(pj / K*j) has to decrease by at 
least a factor 2 for it to be competitive with the constraints from AM^ and AMb„- This requires a 
dedicated effort from the Lattice community, which is already under way HS^I We hope that a robust 
calculation of SU(3) symmetry breaking in £ will soon be undertaken, as experiments are fast approaching 
the SM-sensitivity in R(pj/K*j). 

6.2 A Global fit of the CKM parameters and predictions for a, 7 and AMb s 

To conclude this section, we give here the allowed ranges of the CKM-Wolfenstein parameters and the 
angles of the unitarity triangle obtained from a global fit of the data. Several input quantities that enter 
in the fits have evolved with time and their current values differ (see, Tabic |3J from the ones given in 
the CERN CKM- Workshop proceedings^, and also from those used in the popular fits of the CKMfitter 
group 1221. Hence, a consistent update is not out of place. 

First, a few words about the fits. We follow the Bayesian analysis method in fitting the data. 
However, it should be pointed out that the debate on the Bayesian vs. Non-Bayesian interpretation of 
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Table 3. Input parameters used in the CKM-unitarity fits. Values of the other parameters are taken from the recent 
PDG review. El 



Parameter Input Value 

A 0.2224 ± 0.002 (fixed) 

|Vci,| (41.2 ± 2.1) x 10~ 3 

1Kb | (3.90 ± 0.55) x 10- 3 

a^K s 0.736 ± 0.049 

\e K \ (2.280 ± 0.13) x 10- 3 

AM Bd 0.503 ± 0.006 (ps) -1 

r)i(m c (m c ) = 1.30 GeV) 1.32 ± 0.32 

r)2 0.57 ±0.01 

% 0.46 ±0.05 

m c (m c ) 1.25 ±0.10 GeV 

m t (m t ) 167 ±5 GeV 

B K 0.86 ±0.15 



f Bd y/B^ (210±24)MeV 
T] B 0.55 ±0.01 

£ 1.19 ±0.09 

AM Bs < 14.4 (ps) _1 at 95% C.L. 



data is a lively subject and it has implications for the CKM fits. In the present context the issues 
involved and the quantitative differences in the resulting profiles of the unitarity triangle are discussed 
in depth in the literature | 60 | 9 | 214 | These differences are not crucial for our discussion, but the input 
values of the parameters are indeed important. To incorporate the constraint fro m AM g. , we have used 
the modified-;^ 2 method (as described in the CERN CKM Workshop proc eedings B^l) , which makes use 
of the " Amplitude Technique" introduced by Moser and Roussarie I I in which the time-dependent 
oscillation probabilities are modified to have the dependence P(B°(0) — > B®(i)) oc (1 + Acos AM Bs t) 
and P(B°(0) -> ~B°(i)) oc (1-^4 cos AM B t). The contribution to \ 2 of the fit from the analysis of AM Bs 
is obtained using the following expression 



,Y 2 =2 



(150) 



where A and a a are the wor ld average amplitude and the error, respectively. Relegating the details to a 
forthcoming publication 1 ^ 1 , the main results are summarized below. 

The constraints in the (p,fj) plane resulting from the five individual input quantities (Rb, £k, AM Bd , 
AM Bs , and a^K s ) are shown in Fig. 1111 and correspond to 95% C.L., in contrast to the ones shown 
in Fig. ^3 The resulting 95% C.L. fit contour is shown in Fig. The apex of the triangle for the 
best-fit solution (x 2 = 0.57 for two variables) is shown by the black dot and corresponds to the values 
(p,fj) = (0.17,0.36), with the 68% C.L. range being 

0.11 <p< 0.23, 0.32 < < 0.40 . (151) 
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We also show the constraint from the 90% C.L. upper bound on R(frf / 'K*j) for the value C = 0.75, 
though we have not used this input in the fits. The 68% C.L. ranges of the Wolfenstein parameters 
A, p and fj, together with the corresponding ranges for the CP-violating phases a, (3, 7, and the mass 
difference AA1b s are given in Table 0| Note that the fit-range for sin 2/3 coincides practically with the 
experimental measurement sin 2/3 = 0.736 ±0.049. If we take away the input value of sin 2(3 from the fits, 
and instead determine sin 2/3 from the unitarity fit, we get sin 2/3 = 0.730 ±0.085, which is in remarkable 
agreement with the experimental value, but less precise. In fact, similar estimates of sin 2/3 from the 
CKM unitarity const raints were obtained in the SM by several groups long before its measurements from 
the CP asymmetries I ^17 | 218 | 219 | 220 | 221 | 222 | Now, that the measurement of a^K s lS quite precise and 
its translation into sin 2/3 is free of hadronic uncertainties, this input has reduced the allowed parameter 
space in (p, ?7)-plane - a feature already noted in the literature EDIIaSEEIll 

Prediction of AMg e from the fits deserves a discussion. First of all, if we take away the bound 
on AMg s given in Table [21 from the fits, the allowed range for p becomes lot larger though fj remains 
essentially unchanged, yielding 0.08 < p < 0.27 and 0.31 < fj < 0.41 at 68% C.L. The corresponding 
allowed range for AMb b in this case is 13.0 < AMb s < 21.6 (ps) _1 , with the 95% C.L. interval being 
[8.6,26.2] (ps) _1 . This range is to be compared with the corresponding one [10.2,31.4] (ps) _1 from the 
CKMfitter group ^01. The reason for this apparent mismatch lies in the values of the input parameters £ 
and fB d \/BB d , for which we take the currently updated Lattice values 1.19 ± 0.09 and (210 ± 24) MeV, 
as opposed to the values 1.21 ± 0.04 ± 0.05 and (228 ±30 ±10) MeV, respectively, used by the CKMfitter 
group. With their input values, the resulting 95% C.L. range for AMb s from our fits becomes [10.7, 30.3] 
(ps) -1 , which is almost identical to theirs. Hence, it is important to know £ and fs d \J Bs d more precisely. 
Including the AM Bs bound in the fits, we get at 68% C.L. 14.4 < AM B „ < 20.5 (ps)" 1 with the 95% 
C.L. range being [14.4,23.7] (ps) _1 . 

Table 4. 68% C.L. ranges for the Wolfenstein parameters, CP-violating phases and AMg s from the CKM-unitarity fits. 



Parameter 68% C.L. Range 



p 


0.11 - 0.23 


V 


0.32 - 0.40 


A 


0.79 - 0.86 


sin 2/3 


0.68 - 0.78 


P 


21.6° - 25.8° 


sin 2a 


-0.44 - 0.30 


a 


81° - 103° 


sin 27 


0.49 - 0.95 


7 


53° - 75° 


AM Bs 


14.4 - 20.5 (ps) 



Concluding this section, we remark that the overall picture that has emerged from the current 
knowledge of the CKM parameters and hadronic quantities is that the CKM theory describes all data 
on flavour physics remarkably consistently, including \ex\ and a^,K s - Hence, it is very likely that CP 
violation in hadronic physics is dominated by the Kobayashi-Maskawa phase. Despite this impressive 
synthesis of flavour physics, a clean bill of health for the CKM theory still awaits a number of benchmark 
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Figure 11. Constraints on the p - fj plane from the five measurements as indicated. Note that the curves labelled as 
R(p-f / K*f) and AMg g are obtained from their 95% C.L. upper limits 0.047 and 14.4 (ps) -1 , respectively. The fit contour 
correspon ds to 95% C.L., with the fitting procedure explained in the text. The dot shows the best-fit value. (Figure updated 
from Ref.EHU) 



measurements. These include quantitative determination of the other two angles a (or fa) and 7 (or fa), 
and AMb s . In all these cases, experiments have well-defined targets to shoot at, as the steady progress 
in the knowledge of the CKM parameters, and quite importantly the precise measurement of sin 2(3, have 
reduced the allowed parameter space substantially. The 95% C.L. ranges for these quantities resulting 
from the fits described earlier are: 

70° < a < 114° , 43° < 7 < 86° , 14.4 (ps)" 1 < AM Bs < 23.7 (ps)" 1 . (152) 

The long run to these goal posts has already started. We anticipate significant measurements of all 
three within the next couple of years in experiments at the B factories and Tevatron, but definitely in 
experiments planned at the hadron colliders (LHC-B, ATLAS, CMS and BTEV), which will measure all 
three angles (a,/?, 7) accurately, as well as AMb s and a number of rare -B-decays, such as B s — > /i + /J~. 
Present situation together with some theoretical suggestions is reviewed in the next section. 

7 CP Violation in B-Meson Decays 

In the preceding section, we briefly discussed the CP asymmetries \en\, Re (e'/e) and a^K S i representing 
three different ways in which CP violation has been mea sured so far in the weak decays of the hadrons. 
However, CP violation is a rich and diverse phenomenon E23I This is illustrated in Fig.^J showing its 
various manifestations in the decays of a neutral meson (P°) into a final CP eigenstate (fcp)- Here, 
P° stands for any of the mesons K°,D°,B^,B^. As CP asymmetries arise due to the interference of 
two different amplitudes with their own weak phases, this figure shows that there are three generic 
manifestations of CP violation in P° decays. 

1. Interference between the two decay amplitudes, called A\ and A^. 
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Figure 12. CP violation in neutral mesons decaying into a final CP eigenstate. Each directed line represents an amplitude 
with its own weak phase, and the connecting double arrows between these lines indicate possible interference patterns 
involving these amplitudes. (Figure taken from NirESl) 



2. Interference between the mass (M 12 ) and width (r 12 ) parts of the P° - P° mixing amplitude, giving 
rise to a non- vanishing relative weak phase. 

3. Interference of the decay amplitudes with and without mixing, P° — > P° — > fcp and P° — > fcp. 

In the decays of charged mesons P^ (such as B ± , D ± , Df ) and baryons, CP violation can take place 
only through the interference in the decay amplitudes. 

Assuming CPT invariance, time reversal violation, or T-violation, implies CP violation. In the KM 
theory, CP- and T-violation have a common origin, namely the KM-phase of the CKM matrix. T-violation 
has been established in the K°-K° system. The measured T-violating asymmetry At = (6.6±1.0) x 10~ 3 
is consistent with the measurements of the CP- violating parameter ex ™, hence with the KM phase. T- 
violations have also been searched for in flavour conserving transitions. The current b est u pper limits 
on the T -violating electri c dip ole moments (EDMs) have been obtained for the neutron 12251 Thallium- 
205^ and Mercury-199^1 

\d n \ < 6.3 x 1CT 26 c cm (90% C.L.) , 

K(T1- 205)| < 1.6 x 1CT 27 e cm (95% C.L.) , (153) 
|4(Hg - 199)| < 2.1 x 10~ 28 e cm (95% C.L.) , 

where the last two are to be interpreted as limits on the EDM of the electro n. Ju dging from the current 
upper limit on d n and the prediction of the KM theory d n = O(10~ 33 ) e cm EMI T-violations in flavour 
conserving transitions do not provide any information on the KM phase. As a measurement of d e is not 
foreseen in the SM either, a positive result in any of the three EDMs will be a proof that additional weak 
phases are operative in the weak interactions of the hadrons and leptons. In particular, Supersymmetry 
has a myriad of weak phases, including the phases of the A and /x terms, which in some parts of the 
supersymmetric parameter space are large enough to yield values of the EDMs of the neutron, Hg-199 
and Tl-205 atoms just below the present experimental upper bound E2HI J n some cases, supersymmetric 
weak phases will also manifest themselves in CP violation in flavour changing transitions in the B- and 
JC-meson secto rs. For example, such theories may lead to CP-viola ting effects in the radiative decay 
B — > X s 7 E2QI Recent analyses incorporating t he co nstraints from rfn, E2U predict CP asymmetry in this 
decay close to the current experimental bounds 12221 

-4cp(A s 7) = 0.004 ± 0.051 ± 0.038 => -0.107 < Acp(X s j) < 0.099 (90% C.L.) , (154) 
A C p(K*-f) = (-0.5 ± 3.7)% -0.066 < A C p(K*-f) < 0.056 (90% C.L.) . (155) 



43 



The bounds on Acp(K*"f) are stronger but they still allow this CP asymmetry to be of 0(5%), much 
larger than the SM-based expectations l ^33 | 2 34| |./4cp(X s 7)| < 0.5%. It is entirely conceivable that 
precision studies of CP violation in B- and i-T-decays may require the intervention of some of the weak 
phases in supersymmctric theories. The search for non-KM weak phases is an important part of the 
research programme at the i?-factories and later at the hadron colliders and has to be pursued vigorously. 
However, we will not discuss these scenarios here, as the principal focus of this review is on the CKM 
phenomenology and current data show no significant deviations from the CKM theory. 

In this section, we first discuss each of the three classes of CP asymmetries depicted in Fig. ^| and 
define the underlying physical quantities which are being sensitively probed in each one of them. This 
will be followed by a summary of the existing results on CP asymmetries from BABAR and BELLE in 
some of the main decay modes of interest for our discussion, such as J /ipKs, 4>Kg 1 rj ' K$, tttt, Kit. In the 
SM, these asymmetries provide information on the weak phases a, (3, 7 (or <^i, </>2, ^3), and we review their 
current knowledge from the £>-factory experiments. 



7.1 CP violation in decay amplitudes 

As opposed to if-mesons, direct CP violation in the £>-meson sector is pote ntial ly a very prolific phe- 
nomenon as the number of decay channels available in the latter is enormous l ^*H j n p rac tice, how ever, 
measuring some of these asymmetries, which are estimated at several percent or somewhat higher level! 
in decay modes with typical branching ratios of order 10 -5 , requires O(10 9 ) BB pairs (or more). With 
the present integrated luminosity of O(10 8 ) BB mesons recorded by the BABAR and BELLE detectors, 
there are now indications in the data that these experiments are obs erving the first direct CP asymmetry 
in B decays. The case in point is _4cp(A" + 7t~) = (—9.5 ± 2.9)% E^l having a 3<r significance. 

To study direct CP asymmetry, we note that unitarity can be used to write the decay amplitudes 
A(B — > /) and its CP-conjugate A(B — > /) as: 



A(B -> /) = \A 1 \e +Wl e +tSl + \A 2 \e +t02 e +lS2 , 

(156) 

A(B -> /) = \A 1 \e- t9l c +tSl + \A 2 \e- W2 e +U2 , 

where 9\ and 9 2 are the weak phases, and the strong interaction amplitudes |Ai|e +!l5i (i — 1,2) contain 
the CP-conserving strong phases 5j. With the help of these amplitudes, the CP-rate asymmetry can be 
written as 

T(B^ f)-T(B^ f) _ 2|A 1 p 2 |sin(<y 1 -fe)sin(fl 1 -0 2 ) 
*CPV)- T{B ^ f) + T(B^ f) \A 1 \ 2 + 2\A 1 \\A 2 \cos(S 1 -6 2 )co S (e 1 -e 2 ) + \A 2 \^ { 1 

The goal is to extract the weak phase difference 9 = 9\ — 9 2 from the measured partial rate asymmetries, 
for which knowledge of the strong interaction amplitudes \Ai\, \A 2 \ and the strong-phase difference 
S = 5± — S 2 is essential. The required strong phase difference involving a tree and penguin amplitudes 
5 = Sp — St, or involving two penguin amplitudes with different strong phases 5 = Sp — 5p> , are generated 
by the so-called Bander-Silverman-Soni mechanismESZI i n addition, final state interactions also generate 
strong phases, which have to be estimated or measured. 

From the CP-averaged decay rate (T(B — > f) + T(B — > f))/2 and the CP asymmetry „4cp (/), one has 
two equations with four unknowns which we take as r = |^4_2 1 / |^4i U an d I-^iIj with r < 1. Hence, one 
has to use additional experimental input and assumptions. For the decays B — > h\h 2l with hi belonging 
to a E/(3) nonet of (usually vector or pseudoscalar) mesons, an approximate SU(3) symmetry is invoked 
which allows to express |j4i| in terms of a known branching ratio dominated by this amplitude. This 
then leaves only three unknowns with two constraints, which is still not sufficient to determine all the 
parameters, but allows to derive bounds on the weak phase 9 which could be useful if data is benevolent. 

While the phenomenology of direct CP violation is rich USED a theoretically robust description is 
difficult and not yet at hand. Determining the weak phases from the observed asymmetries and decay 
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rates requires a good control over the ratios of the amplitudes \Ap/At\, \Ap>/At\, and \Ap/Ap/\ etc., 
and the corresponding strong-phase differences 5p — St, §p> — St and Sp — 8 pi, where we are assuming 
that there are two different penguin amplitudes, one generated by the exchange of gluons and the other 
by electroweak bosons (7 and Z). Of course, direct CP asymmetries also arise when two different tree 
amplitudes interfere, such as in Acp(DK). Anyway, calculating them from first principles with a non- 
perturbative technique, such as Lattice QCD, is simply not on the cards, as the amplitudes |^4i|e I<5i 
depend on (hxh2\0\B), where O is a four-quark operator. We know all too well that the vengeance 
of strong interactions is merciless in the analysis of Re(e'/e) in if-meson decays involving the matrix 
elements (int\0\K l s) 1 which has hindered the extraction of any information on the CKM parameters 
from this measurement. 

One would like to argue that strong interaction effects are tractable in B-decays due to the large 
mass of the &- quark, and hence calculable using techniques based on 1/ mf,-expansion and perturbative 
QCD 05212221 However, the dynamical scale in B — > h\}i2 decays is not set by the inverse &-quark mass, 
1/mb, but by a lower value, of order 1 -2 GeV - the typical virtuality of a gluon nested somewhere in 
the Feynman diagrams. This raises the question if perturbative methods are adequate at such low scales. 
Jury is still out on this issue, but when the jury returns one should not be surprised at an unfavourable 
verdict for the perturbative-QCD inspired factorization models. 

Some help in extracting the weak phases from data can certainly be sought by using arguments based 
on the isospin and flavour SU(3) symmetries. The most ce lebra ted use of the isospin symmetry is in the 
extraction of the phase a (or ^2) from the B — ► mr decays 12221 While isospin is a good symmetry, there 
are still missing experimental links to complete the chain of arguments to determine the weak phase a in 
a model-independent way. However, use of the SU(3) symmetry is more problematic. We have seen that 
the issue of SU(3)-breaking in simpler quantities, such as the ratios £ and f, is far from being settled. 
In the much more complicated situation in non-leptonic decays, SU(3) symmetry breaking effects are at 
best modeled, often based on the assumed properties of the factorized amplitudes. This lack of a robust 
theoretical basis and/or data introduces hadronic uncertainties in the extraction of the weak phases from 
data. We shall discuss applications of these methods in the context of B — ► tttt and B — > Kir decays to 
quantify some of the issues involved. 

7.2 Indirect CP violation involving B° - B° mixing amplitudes 

Indirect CP asymmetries involve an interference between the absorptive and dispersive parts of the 
amplitudes in the B° d - B° d and B° - ~B° mixings. Their experimental measures are the charge asymmetries 
-4 S i(-B°) and _4 sl (B,°): 

V ; T(B°^£+X)+T(B°^£-X) 1 + \q/p\ 4 



where the ratio q/p is defined as follows 



g\ 2 _ 2M* 2 -iT* 12 
2Af 12 - iT-ii 



(159) 



and the off-diagonal elements M12 and Ti2 govern the mass difference (AMp) and the width-difference 
(ATp) between the two mass eigenstates, respectively. Thus, 

Parameterizing the ratio Ti2(B q )/Mi2(B q ) = r q e lt * q , SM estimates are 

A SL (B° d ) =r d sm( d ~O(10- 3 ), A SL (B° S ) = r s S in( s ~ 0(10^) . (161) 
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Present measurements yield ESI „4 sl (jB°) = (0.1 ± 1.4) x 10~ 2 . T he bound on Asl(B^) does not probe 
the SM, but constrains some beyond-the-SM (BSM) scenarios l ^4U | Combini ng th e data on the direct CP 
asymmetry Acp(B + — > J/ipK + ), one can get an improved bound on \q/p\ 12231 anc [ we shall discuss it 
later. However, for all practical purposes discussed here, one can set \q/p\ = 1 for the -B^ system. 
Currently, there is no experimental bound on Asl(B®). 

It is worth pointing out that the analogue of the charge asymmetry 1158(1 in the K° - K° system, 
namely the CP-violating asymmetry in the semilcptonic decays Kl — > Tr^l^Vg, defined as, 

_ T(K L -> n-e+v t ) - T(K L -> tt+Tj^) 
£ ~ T{K L -» tt-£+^) + L(X L -» tt+^-z/^) ' 1 j 

is a well measured quantity 121 § e = (3.27 ± 0.12) x 10 -3 . This value is consistent with the relation 
5( = 2Re(ei<-)/(l + \ck\ 2 ) — 2Re(e^), and the experimental value^of Re(e/f). Hence, the asymmetry 6e 
arises entirely from the imaginary part of M\2(K), which in the SM is given by lm(V t *Vtd)- 

Current and planned experiments at the e + e~ and hadronic B factories are not anticipated to reach 
the SM-sensitivity in Asl(B^) (even less so in Asl(Bs))> an d ; hence, a measurement of any of these 
asymmetries will be a sure signal of BSM physics. 

7.3 Interplay of mixing and decays of B° and B° to CP eigenstates 

As already mentioned, this class of CP asymmetries involves an interference between the decays B° — > fcp 
and B° — > B° — > fcp- Due to mixing these CP asymmetries are time-dependent, and the two time- 
dependent functions cos(AAfei) and sin(AMst) can be measured, allowing to extract their coefficients 
Cf and Sf, respectively: 

A tC p(t) = r W f )- > /)~ r ( g0 ( f )- > /) = _ c C0S (AM B t) + S f sin(AM B t) , (163) 

where 

Cf ~ ~ Af - iTJXff Sf -iTJXfT' (164) 

and the dynamical quantity Xf is given by: 

\f = (q/p)p(f); M = ^f)- ( 165 ) 

Concentrating now on B® and B^ decays, the amplitudes A(f), A(f) and q/p are defined as follows^ 

A{f) = ; !(/) = ; q/p = ^ = e ~ 2 ^™ = e" 2 ^ . (166) 

If only a single decay amplitude is dominant, then one can write: 

p(/)=ry/e- 2i ^ =► |p(/)|=l, (167) 

where rjf = ±1 is the intrinsic CP-Parity of the state /. In this case, |A/| = 1, and one has 

S f = lm(Xf) = -r)f sin2((/> mixing + </> d ccay) , C/=0, (168) 

and the CP asymmetry .Afcp has a very simple interpretation: 

Atcp(t) = S f sm(AM B t) . (169) 

If, in addition, </>dccay = 0, which is the case for the transitions b — > ccs, b — > sss and b — ► dds, then a 
measurement of Sf is a measurement of sin2# m i X i ng (modulo the sign —rjf). The current thrust of the 
BABAR and BELLE experiments in CP asymmetry measurements is indeed in extracting Sf and Cf for 
a large number of final states. 
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7.4 Current status of CP asymmetries in b — > ccs and b — > sss decays 

We now discuss the coefficients S/ and C/ for / = J/ipKs,4>Ks,r]'Ks,K + K~~K s . Current data on 
these and some other final states are summarized in Table [5] The values quoted for f = J/ ipKg are 
averages over several decay modes of the quark decay b — ► ccs and include also the final states J/ ipK L 
and J/ipK*, with the latter angular momentum analyzed into states wi th w ell-defined CP-parities and 
taking into account the intrinsic CP-parity and some other related states^^^. Also, the state K + K~ Ks 
is not a CP eigenstate and an angular analysis like the one carried out for the B — > J/ ipK* case has not 
yet been undertaken due to limited stati stics. However, using arguments based on branching ratios of 
related modes, BELLE collaboration IMH concludes that the K + K~K$ state is predominantly a CP-even 
eigenstate, with a fraction 1.03 ± 0.15 ± 0.05. 

In the SM, Sj/^k s an d Cj/j/>if s are determined by the tree amplitude b — > ccs with the penguin 
amplitude suppressed by A 2 . In any case, both the T and P amplitudes have the same phase, and hence 
in the SM Sj/^Ks = sin2/3 (or sin20i) and Cj/^k s = ( or \^J/-ipK s \ — !)■ Comparison of Sj/^k s with 
the indirect estimates of the same discussed earlier shows that the agreement with the SM expectations 
in this decay mode holds quantitatively. This is a great triumph for the KM mechanism of CP violation. 

The present measurement Cj/^,xs = 0-052^q'q46 ^ s m agreement with no direct CP violation in the 
decays i? ± — > J/ipK^. The current bound 

Acp(B+ -> J/^K+) = - -0-007 ± 0.019 , (170) 

yields \A/A\ = 0.993 ±0.018. Combined with \q/p\ = 0.9996±g;gg|f from A SL (B°), this yields \\j/^k s \ = 
0.992 ± 0.019, in precise agreement with \Xj/^k s \ = !• 

Going down the entries in Table[3 we note that the final states (/)Kg and K + K~Kg are determined by 
the penguin transition b — > sss, which has the weak phase n. The final state rj' K$ receives contributions 
both from the penguin b — > sss and b — > sdd amplitudes and the tree amplitude b — > uus, due to the 
uu content of the 7/-meson wave function. The tree amplitude, however, is CKM suppressed. So, to a 
good approximation, also this final state is dominated by the penguin transitions b ~> sss and b — > sdd, 
and has the weak phase it in the decay amplitude. Thus, in the SM, we expect that for these decays 
(f = J/^K s ,4>Ks,v'Ks,K+K-K s ) 

-VfSf = sin 2/3, C f = 0. (171) 

The curr ent m easurements of Sf a nd C f are summarized in Table [S] The data are from BELLE l^iil and 
BABArEEJ updated by Browder^ni 

at the Lepton-Photon Conference and the Summer 2003 updates 
by HFAG'^I Note that ru^K s — V^Ks = Vri'K s = — 1, but rj K + K - Ks = +1 due to the dominance of the 
CP = +1 eigenstate in this mode. One sees that within the experimental errors, SM predictions Cf = 
for the final states specified in (|171|1 are in agreement with the data, though the errors in the individual 
modes are still quite large. On the other hand, S^k s an d to a lesser extent also S v >k S i appear to be out 
of line with the SM expectations l|171l) . Howeve r, it shou ld be noted that the two experiments B AB AR 
and BELLE are not consistent with each other I ^P I ^ w j t h s^ Kg (BELLE) = -0.96 ± 0.50±£S? and 
S^Ks (BABAR) = +0.45 ± 0.43 ± 0.07, differing by 2.1a. Following the PDG rules ^, one has to scale 
the error in S^Ks given in Table|S]by this factor, yielding S^Ks = — 0.14 ±0.69, which is only 1.3cr away 
from Sj/^pxs = 0.736 ± 0.049, and hence the difference between the two is not all that compelling. 

As the current statistics is low, it is helpful to combine the three (b — ► sss) p engu in-dominated final 
states. Defining —nfSf = sin2/3 e ff, and averaging over the three final states gives'^ 

(sin 2/3 cff ) = 0.24 ±0.15 (C.L. = 0.11) , 

(172) 

(C off ) = 0.07 ± 0.09 (C.L. =0.76). 

The two experiments (BABAR and BELLE) are still hardly compatible with each other in (sin2/3 c ff); 
in contrast the measurements of (C e g) are quite consistent. The current average of (C e g) within errors 
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Table 5. The coefficients Sf and Cf = —A-f from time-dependent CP asymmetry A{cp(t), taken from the HFAG listingsE^l 
[BE(BA) stands for a BELLE (BABAR) entry]. 



c f 



c 



J/4>K S 



= 0.052 



+0.048 
-0.046 



C^ Ks = -0.04 ±0.24 
A' s = 0.04 ±0.13 
Ck+k-k, [BE] = 0.17 ± 0.16 ± 0.05 



C 



= 0.13 ±0.24 



C D ^[BA] = 0.28 ± 0.23 ± 0.02 
C + -(D*+D-)[BA] = -0.47 ±0.40 ±0.12 
C-+(D*- L>+)[BA] = -0.22 ±0.37 ±0.10 



-0.38 ±0.16 



Sr 



Sj/^k s = 0.736 ± 0.049 



Scf>K s = 
S V 'K S ~ 
Sr+K-Ks I be ] 



-0.14 ±0.33 
+0.27±0.21 
= -0.51 ±0.26 ±0.05 



Sj/^o = -0.40 ± 0.33 
S D , W [BA] = 0.06 ± 0.37 ± 0.13 
S+-(D*+D- )[BA] = -0.82 ±0.75 ±0.14 
S- + (D*-D+)[BA] = -0.24 ±0.69 ±0.12 



-0.58 ±0.20 



also agrees with the corresponding quantity measured in the b — ► ccs transitions, with Cj/^k s ~ (Ceff) = 
0.018 ±0.10. However, the current measurements yield 121121521 

sin 2/3 - (sin 2/3 cff ) = 0.50 ± 0.16 , (173) 

which is a 3.1er effect on the face value. Taking into account the scale factor in (sin2/3 c ff) increases the 
error, yielding sin 2/3 — (s'm2(3 c g) — 0.50 ± 0.25, which differs from by 2a. This difference is more 
significant than Sj/^Ks ~ S<$>K S i but still does not have the statistical weight to usher us into a new era 
of CP violation. We have to wait for more data. 

While there are already quite a few suggestions in the recent literature explaining the difference 
S^k s — Sj/^Ks m terms of physics beyond the SM, we will not discuss them as the experimental sig- 
nificance of the effect is marginal. However, a more pertinent question to ask is: How well are th e SM 
equalities given in (|171fl satisfied? This point has been investigated recently by Grossman et al. ESI 
Following their notation, the SM amplitudes in these decays can be parametrized as: 

A f = A(B° -> /) = V; b V cs a c f + V: b V us a u f , (174) 

where is dominated by the b — > sss gluonic penguin diagrams and gets contributions from both 
penguin and b — > uus tree diagrams. The second term is CKM-suppressed compared to the first: 



Im 



V*V 



v: h v us 



V* V 



sin 7 = 0(A 2 )sin7. (175) 



It is conceivable though not very likely that the CKM-suppression is offset by a dynamical enhancement 
of the ratio a'j/a'j. Note that \a u f / a c f \ ~ 1 (from penguins), but \aJ/a C f\ (from tree) could be ^> 1. To 

quantify this, Grossman et al. ' ' define 

V* V a u 
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SU(3) allows to put bounds on 

— rjfSf — sin2/3 = 2cos2/3sin7COS<5/|£/| , 

(177) 

Cf = — 2 sin7sin , 

where 6f = axg(aj / a*) and £f also characterizes the size of Cf. Note that 6f can be determined from 
tan (5/ = (rjfSf + sin 2/3) /(Cf cos 2/3). However, present bounds are not very restrictive due to lack of 
information on so me decays and additional assumptions are required to be more quantitative. Typical 
estimates are EI31 

\£r,'K s \ < 0.36 [SU(3)]; 0.09 [SU(3) + Leading N c ] , 

\^k s I < 0.25 [SU(3) + Non - cancellation] , (178) 
\tK+K-K a \ <0.13 [C/-Spin], 

where the various assumptions in arriving at the numerical inequalities have been specified. More data 
and further theoretical analysis are required to further restrict in a model-independent way. Hence, 
at this stage, one should use values of indicated in (|178|) . 

7. 5 Current status of the CP asymmetries in b — > ccd decays 

CP asymmetries in some of the final states which are induced by the transition b — > ccd have 
also been studied by the BABAR and BELLE collaborations. The coefficients Sf and Cf for / = 
J/ijjir , D*D*, D* + D~, D*~D + measured through the time-dependent CP-asymmetry in these chan- 
nels are given in the lower part of Table [SJ Note that the entry for / = J/tpir is the average of the 
BABAR EH and BELLE experi ments b ut the entries for / = D*+D-, f = D*~ D+ and / = D*D* 
are from the BABAR collaboration M 4 o | ^ 4 I \ gj one w ith the BELLE results not yet available. Just like 
the b — > ccs case, the tree amplitude in the b — ► ccd transition does not have a weak phase. However, 
as opposed to the b — > ccs case, where the penguin amplitude has the weak phase n, now there could 
be a non-negligible contribution from the b — > d penguin amplitude, which carries the weak phase (3. In 
terms of the CKM factors, both the T and P amplitudes are of order A 3 . The tree-penguin interference 
with different weak phases may lead to direct CP violation, giving rise to Ct ^ 0, or |A/| ^ 1. Also, 
in this case, the equality —Sj/^o — Sj/^k s or — S D «jp-(+) — Sj/^k s e ^ c - can be violated. Here the 
sign(+) stands for the CP= +1 component of the Vector- Vector final state, which for D*D* dominates 
the decay with the CP-odd component given by^^R± = D*IF(P w&ve)/D*W = 0.063±0.055±0.009. 
Present data are consistent with the SM and no direct CP violation in these decays is observed. Also, the 
measured coefficients 5*/ for these decays do not show significant deviations from the SM expectations 
(ignoring penguins) —rjfSf — Sj/^^s- This brings us to the entry in the last row in Tablc[S]for / = 7T7t, 
which we discuss below together with the two methods most discussed in the literature to determine the 
angle a (or fa). 

7.6 Measurements of the weak phase a (or fa) in B-meson decays 

In the previous section we have given the 95% C.L. range for the weak phase a obtained from the CKM 
unitarity fits: 70° < a < 115°. This phase will be measured through CP violation in the B — * irir and 
B — > (07T decays. To eliminate the hadronic uncertainties in the determination of a, an isospin analysis 
of these final states (as well as an angular analysis in the pn case) will be necessary. However, at a less 
rigorous level, data from B — > Kit decays may be combined with the available data on the B — > nn 
decays to extract a from the current data. We discuss both of these methods below. 

The decay B — > 7r + 7r - involves tree and penguin contributions with different strong and weak 
phases. Denoting the strong phase difference by d = Sp — St, the amplitudes can be written as 

- A(B° -> 7T+7T-) = \T\e l ~< + \P\e lS , 
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_ (179) 
-A(B° -> tt+tT) = \T\e~^ + \P\e lS , 



where the T and P components have the CKM dependence (using the Gronau-Rosner convention 
given by V* b V u d and V* b V c d, respectively. The time-dependent CP asymmetry is given by the expression 

AtttC*) = -C™ cos{AM Bd t) + sm{AM Bd t) , (180) 
and the coefficients and arc defined as in l|164|) . with 

A = „ e-^)i±ra^ (181) 

Thus, apart from f3 (or </>i), which is now well measured, we have three more variables, the ratio \P/T\, 
5 and 7 (or 03 ). As this expression stands, it gives information on 7 and not on al However, if the 
penguin contribution were absent (or small), then using the relation a + [3 + 7 = 7r and rj^ = +1, one 
has A^tt = e~ 2l (P +J ) — e 2ra , yielding = and = sin 2a. Hence, the folklore: SVtt measures 
sin 2a. Now, with strong hints from data that \P/T\ is si gnific ant (for example, T-dominance would 
have 2r(7r+7r°)/r(7r+7r-) = 1, but the latest BELLE data 1^1 gives 2.10 ± 0.58 ± 0.25 for this ratio) 
one interprets S„ n as S^r = \/i — sin 2a c s , where both CV^ and a c g involve non-trivial hadronic 
physics. Hence, a transcription of a e ff into a (or 7) is not easy to accomplish in a model-independent way. 

To get a model- independent determination of a (or ^2), one has to carry out the isospin analysis of 
the B — > 7T7T decays suggested by Gronau and London . Defining the various amplitudes as follows: 

A+- = A(B Q -» 7T+7T-); A+- = ^ Tr+O , 

A 00 ee A(B° -» ttV ); A 00 ee A(W -» tt (182) 
A +0 ee A(B+ tt+tt ); ee t^tt ) , 

isospin-symmetry leads to the following triangular relations l ^^ l 

1 A +- + A oo = A +o J_ A+ - + A m = A -o 

V2 V2 

\A+°\ = \A-°\. (183) 

Here, the last equality results from the observation that the amplitudes A +0 and A~° describe decays into 
pure isospin-2 states and do not receive contributions from the QCD penguins, and electroweak pen guins 
may be ignored, as their contribution in the nn system is not expected to exceed a few percent ' '. 
One can also include the contribution of the electroweak penguins in this analysis by using isospin 
symmetry ^ ^, which relates their contribution to the tree amplitudes and eliminate any residual hadronic 
uncertainty in the determination of the weak phase a. 

The two triangles written in the first line in (|183|l have the same base (due to the second line in 
l|183p l and the mismatch in the apex of the two triangles then determines the difference 29 = 2(a c ff — a). 
The determination of 2a goes along the following lines: From the relative phase of the amplitudes A H 
and A^ one gets 2a c g . From the relative orientation of the amplitudes A +0 and A^ one gets an angle 
and finally from the relative orientation of the amplitudes A~° and A^ one gets an angle $. The 
angle 2a is then obtained from the difference 2a = 2a c g — $ — $. 

From the magnitudes | A^ | and \A^ | of the six amplitudes given in the isospin relations l|183|) , the only 
missing pieces in the experiments are |A 00 | and |A 00 |. However, through the measurement of the charge- 
conjugate averaged branching ratio B a /B° — > 7r°7r°, the combination |A 00 | 2 + |A 00 | 2 is now known. This 
branching ratio together with some of the other B — > irir and B — > Ktt branching ratios is given in Table 
El where the entries are from the Lepton-Photon 2003 conference review by Fry I I. The Gronau-London 
isospin analysis can not be carried out for the time being. Howe ver, theoretical bounds on 9 (or a g) 
have been proposed. For example, the Grossman- Quinn boundEEHl sin (a — a c g )| < \/r(7r 7r )/r(7r + 7r ), 
with the branching ratios given in Table [S] yields \a — a c g\ < 48° at 90% G.L. I ^" 4 !, and hence currently 
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not very help ful. For other suggestions, see recent papers by London, Sinha and Sinha l zoo l and by 
Gronau et al. ESEl 



Table 6. Summary of branching fractions (xlO 6 ) and Acpfoih?) (Source: Lepton-Photon 2003 reviewESZ 



Decay Mode 


B{hxh 2 ) 


A c Ahih 2 ) (%) 


K+7T- 


18.2 ±0.8 


-9.5 ±2.9 


K°n+ 


21.8 ± 1.4 


-1.6 ±5.7 


K+n° 


12.8 ± 1.1 


0±7 


K n° 


11.9 ± 1.4 


3 ±37 




4.6 ±0.4 




7T+7T 


5.3 ±0.8 


-7 ±14 


7r°7r° 


1.90 ±0.47 





We now discuss the information that can be obtained on the phase a by invoking SU(3) relations 
between the B — > tttt tree an d pen guin amplitudes (T and P) and the corresponding amplitudes in the 
B — ► Kn decays (T' and P') ' I In the present context, one may relate the amplitudes in the decays 
B+ -> K a ir+ and B° -> tt+tt-. Writing 

- A(B+ K\+) = \P'\ e" - \P\ e iS . { K . , (184) 

/tt tan fl c 

where a small term with the weak phase 7 has been neglected and factorization of the decay amplitudes 
is assumed. With the known values of f K , fa, tanfl c and B(B° d -> ir+ir-)/B(B+ -> if°7r+) = 0.23± 0.03 
(from TablcEl, one gets \P/T\ ~ 0.3. This allows to con strain a from the present measurements of 
and Ctttt I2MI ^48 | Another variation on the same theme EHU jg to use the data from the B d — > IsT^^ 
decays instead of B + — * AT°7r + . Again, using flavour SU(3) symmetry and dynamical assumptions, 
one can extract a from the current data on S^^, C nv and B(B d — > 7r + 7r~)/£>(i?° — > if tt^). These 
methods actually give information on 7 (or </> 3 ), as discussed above, and the consequences of the current 



measurements have been recently worked out by Fleischer laSDl getting a value of 7 in the SM-ball park 
with (3 ~ 24°. 

With the average S^tt = —0.58 ±0.20 and C„„ = —0.38 ±0.16, as given in TableEl we see that within 
measurement errors, — SVtt = Sj/^Ks l& n °t violated, which implies no direct CP violation, but 
deviates from by about 2.4 er, which is a signature of direct CP violation. So, at present, the inferences 
from C W7r and S nn are not quite equivocal. The relation — = Sj/^k s IS expected to be violated in the 
SM, as the phases a and (3 are numerically quite different. A large value of \C nn \ would also imply a large 
strong phas e 5, which would put to question the validity of the QCD factorization framework in B — ► irn 
decays ' ^Q \ However, it sho uld b e emphasiz ed th at, just like the coefficient S<pK s , also S nn comes out 
very different in the BABAFtEHland BELLeESI measurements: 5™ (BELLE) = — 1.23 ± 0.41±g;gf and 

SWtBABAR) = -0.40 ± 0.22 ± 0.03, and the average has a C.L. of 0.047^. Scaling the error by the 
PDG scale factor, we get = —0.58 ±0.34 and C W7r = —0.38 ±0.27. With the scaled errors, C nn differs 
from by only 1.4cr. We hope that with almost a factor two more data on tapes, the experimental isuues 
in B — > 7T7T decays will soon be settled. 

Having stated the caveats (theory) and pitfalls (experiments), we use (|184H to illustrate what values 
of a are implied by the current averages S n7r = —0.58 ± 0.20 and C W7r = —0.38 ± 0.16. To take into 
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account t he u ncertainty in the SU(3)-breaking and non-factorizing effects, we vary the ratio \P/T\ in 
the range^^0.2 < \P/T\ < 0.4, which amounts to admitting a ±30% uncertainty on the central value 
of the magnitude obtained in the factorization approach. However, the strong phase S is varied in the 
entire allowed range — ir < S < ir. The results of this analysis are shown in Fig. ^|for four values of a 
which lie in the 95% C.L. allowed range from the unitarity fits: a — 80° (upper left frame), a — 93° 
(upper right frame), corresponding to the nominal central value of the CKM fit, a — 105° (lower left 
frame), and a = 115° (lower right frame), corresponding to the 95% C.L. upper limit on a from our 
fits. We do not show the plot for a — 70°, which is the 95% C.L. lower value of a from the unitarity 
fits, as already the case a = 80° is disfavoured by the current B — > irir data. In each figure, the outer 
circle corresponds to the constraint + — 1; the current average of the BABAR and BELLE 
data satisfies this constraint as shown by the data point with (unsealed) error. The two inner contours 
correspond to the values \P/T\ = 0.2 and \P/T\ = 0.4. The points indicated on these contours represent 
the values of the strong phase 5 which is varied in the interval — ir < S < tt. We see that the B — > tttt 
data is in comfortable agreement with the value of a in the range 90° < a < 115°, with a = 80° outside 
the ±lcr rELIlffG . Also , clS C 7J-7T is negative, current data favours a rather large strong phase, typically 
-90° < S < -30°. This is in qualitative agreement with the pQCD framework 1^1, which yields 
P/T = 0.23±g;g| and -41° < 5 < -32°. (See, also Xiao et al.&B*). 

The case for significant final state rescattering, and hence larg e strong phase differences, in the decays 
B — > PP has also been advo cated by Chua, Hou and Yang motivated in part by the discovery of 
the colour- suppressed decays l 265 l 266 l B° -> D°/i°, where h° = 7r°,?7 or ui, found significantly over 
the factorization-based estimates. An analysis of the current measurements in B — > PP, decays with 
P = tt, K, rj, done with the help of an SU(3) formalism to take into account 8<g>8 — > 8(g) 8 rescattering^H, 
shows marked improvement in the quality of the fit. 

Thus, the study of the CP asymmetry in the B — > tttt decays is potentially the most exciting game in 
town, as precise measurements in this decay mode will not only determine the weak phase a but will also 
decide several important theoretical issues, such as the adequacy or not of the perturbatively generated 
strong phases. To this list of interesting decays, one should add the decays B — > Kir and B — > DK, 
which we discuss below. 



7. 7 Present bounds on the phase 7 from B decays 
The classic method for determining the phase 7 (or (j> 3 ) 

I2fi8l2fi9l270l27ll 

involves the interference of the 

tree amplitudes b — > uW~ — > ucs — > D a K~ and b — > cW~ — > cus — > D°K~. These decay amplitudes can 
interfere if D° and D° decay into a common hadronic final state. Noting that the CP= ±1 eigenstates D± 
are linear combinations of the D° and D° states: Dj. = (D° ± D°)/\/2, both branches lead to the same 
final states B~ — + D±K~ . So, the condition of CP intcrfcromctry is fulfilled. The decays B~ — > D±K~ 
are described by the amplitudes: 

A(B~ -» D±K~) = -= \A(B- D°K-) ± A(B~ -> Wk~)] . (185) 
V2 L i 

Since, the weak phase of the b — * u transition is 7 but the b — > c transition has no phase, a measurement 
of the CP asymmetry through the interference of these two amplitudes yields 7. The four equations that 
will be used to extract 7 are: 

B(B-^D% K~) +B(B+^DlK+) 2 
R± ^ B(B-^DOK-)+B(B+^D°K+) = 1 + ? ' DK ± 2r ° K C ° S &K C ° S 7 ' 



A+ = 



B(B- -> D° ± K-) - B(B+ -> D° ± K+) ±2r DK sinrf D K sin 7 



(186) 



B(B- -» D° ± K-)+B(B+ -> D Q ± K+) 1 + r^ K ± 2r DK cos S DK cos 7 ' 

With three unknowns (^dKj <^dKj 7)1 but four quantities which will be measured R± and A +, o ne has, 
in principle, an over constrained system. Here, tdk is the ratio of the two tree amplitudes I I tdk = 
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Figure 13. Constraints on a from the current measurements Smr = —0.58 ± 0.20 and Cmr = —0.38 ± 0.16 (shown as the 
crossed bars). The outer circle in each frame corresponds to the unitarity condition \G\tt\ 2 + SV-n-l 2 = 1. Theoretically 
estimated range for \P/T\ and the strong phase S varied in the full range — n < 6 < tt are indicated. In the SM, measurements 
of Ctttt and S w7 r have to lie in the region between the two inner contours. This figure shows that values a < 80° are 
disfavoured by the current 7r7r data (upper left frame). The other frames show three allowed values of a within the 95% 
C.L. range from the unitarity fits, which are in agreement with the B — > tttt data. 



\Ti/T 2 \ ~ (0.1 - 0.2), with T x and T 2 being the CKM suppressed (6 -> u) and CKM allowed (b -> c) 
amplitudes, respectively, and <5dk is the relative strong phase between them. The construction of the 
final states involves flavour and CP-tagging of the various D° states, which can be done, for example, 
through the decays D° + — > 7r + 7r~, D°_ — ► A"s7r°, and D° — > K~tt + . Also, more decay modes can be added 
to increase the data sample. 
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Experimentally, the quantities R± are measured through the ratios: 

, , % B(B~ -> D°K-) , , x B(B ± -> Diif±) , 

With all three quantities (R(K/n)andR(Klir)± measured, one can determine R± = 
R(K/ir)±/R(K/ir). Present measurements in the B — > -D-K" and £? — > Z)7r decays have been recently 
summarized by Jawaherey 

ft+ = 1.09 ±0.16, ,4+ = 0.07 ±0.13 [BELLE, BAB AR] , 
R_ = 1.30 ±0.25, A- = -0.19 ±0.18 [BELLE], (188) 
=► r DK = Q-44±°it, (Acp(DK)) = 0.11 ± 0.11 . 

Thus, CP asymmetry in the B — > DK modes is consistent with zero and the ratios R± are both in excess 
of 1. Hence, no useful constraint on 7 can be derived from these data at present. One needs more precise 
measurements of R± and A± to determine 7 from this method. More useful decay modes to construct the 
B — > DK triangle will have to be identified to reduce the statistical errors. In a recent paper, Atwood 
and Soni EZ3 have advocated to also include the decays of the vector states in the analysis, such as 
B~ -> K*-D°, B~ -> K~D*°, and B~ -> K*-D*°, making_use of the D*° -> D°j and D*° -> D°ir° 
modes. This is a promising approach but requires O(10 9 ) BB pairs and reconstruction of many decay 
modes of the D^*'K^ final states to allow a significant measurement of 7. 

A variant of the B — ► DK method of measuring 7 is to use the decays B^ — > DK^ followed by 
multibody decays of the Z?-meson, such as D° — > Kstt~it + , D° —> K $K~ K + and D° — > KsTt~-K + Tt a . 
This was su gge sted some time ago by Atwood, Dunietz and Soni an d was revived more recently by 
Giri et al. ' \ in which a binned Dalitz plot analysis of the decays D° / D° — > Kstt~tt + was proposed. 
Assuming no CP asymmetry in D° decays, the amplitude of the B + — * D°K + — > (Xs 7 '' +7 '' _ )^ + can. be 
written as 

Ai+ = f{m 2 +1 m 2 _) + T DKC ^ +5 ^f{m 2 _,m\) , (189) 

where m 2 , and m?_ are the squared invariant masses of the Kstt + and Kstt~ combinations in the D° 
decay, and / is the complex amplitude of the decay D° — > Kstt + it~ ■ The quantities r^K and 8dk are the 
relative magnitudes and strong phases of the two amplitudes, already discussed earlier. The amplitude 
for the charge conjugate B~ decay is 

M_ = f(m 2 _,m\) + r DK e l{ - 1+5DK) f{rn 2 +1 m 2 _) . (190) 

Once the functional form of / is fixed by a choice of a model for D° — > Kstt + tt~ decay, the Dalitz 
distribution for B + and B~ decays can be fitted simultaneously by the expressions for M + and M_, 
with r^K, Sdk an d 7 ( or ^3) as free parameters. The mo del- dependence could be removed by a binned 
Dalitz distribution Ell 

This method has been used by the BELLE collaboration I^ZsO to measure the angle fa. As the binned 
Dalitz distribution at this stage is limited by statistics and some technical issues involving backgrounds 
and reconstruction efficiencies have to be resolved, an unbinned model-dependent analysis of the B^ — > 
D Q K ± decays followed by the decay D° — > Ksir + Tr~ has been performed. The model for the function 
/ is based on a coherent sum of N two-body plus one non-resonant decay amplitudes, and the N = 7 
resonances used are: K*+ir~, K s p°, K*-ir+, K s lu, K s f (980), K s f Q (U30) and ^(1430)+^-; 

N 

f{m 2 + ,m 2 _) = J2 a ] eUjA j( m l> m -) + be%5 ° > ( 191 ) 

where Aj(m 2 ^ ) m 2 _), a,j and Sj are the matrix element, amplitude and strong phase, respectively, for the 
j-th resonance, and b and Sq are the amplitude and phase for the non-resonant component. Further 
details can be seen in the BELLE paper I I. The result based on 140 fb _1 data has yielded 

r DK = 0.33 ±0.10, fojc = (165115)° , 3 = (92±i?)°. (192) 
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As the errors are quite non-parabolic, they do not represent the accuracy of the measurement. Rather, the 
constraint plots on the pair of parameters (<fo, Sdk) and (tdk, 4>i), which can be seen in the BELLE pub- 
lication^-"^, are used to get the information on these parameters. The resulting 90% C.L. ranges are I t 

0.15 < r DK < 0.50, 104° < 5 DK < 214° , 61° < < 142° , (193) 

and the significance of direct CP violation effect ( including systematics) is 2.4 standard deviations. This 
measurement is in agreement with the SM expectations 43° < 4>z < 86°, though lot less precise. 

As the final topic of this review, we discuss the B — » Ktt decays, which in the SM are dominated 
by QCD penguins. Their branching ratios, averaged over the charge conjugate states, and the present 
measurements of the CP asymmetries are summarized in Tabled The branching ratios listed in this table 
(except for t he B ° — > 7r°7r° mode) are in agreement with the theoretical predictions based on the QCD 
factorization USSl and pQCdGSD approaches. These approaches will be put to more stringent tests with 
precise measurements of the CP asymmetries Acp (Ktt) in various decay modes, and also through the 
ratios of the branching ratios, which provide a better focus on the relative strong phases and magnitudes 
of the QCD penguin, tree and electroweak penguin contributions L££ 6 | 277 | 278 | 279 | 

Concerning the determination of 7 from these decays, whose present status we discuss below, we note 
that the current data on B — > Ktt decays has some puzzling features, which have to be understood before 
we determine 7 from this data. In particular, the following two ratios with their currently measured values: 



R B(K+tt-) =QJ6± 
2B(K°tt°) 

B(K+TT°) 

R c = ^ L = 1.17 ±0.13, 

B{K°TT+) 



(194) 



have received some attention lately l ^80|281| ag p 0SS ibl e harbingers of new physics. It has been argued 
that these data require a much enhanced electroweak penguin contribution than is assumed or inferred 
from the B — > tttt data and approximate SU(3) symmetry. Taking the current data on the face value, a 
range 0.3 < tew = I-PewI/I-^'I — 0-5, and a significant strong phase difference, Sew — St between the 
electroweak penguin and tree amplitudes are n eeded to explain the data. Typical estimates of r^w are, 
however, in the ball park of tew ~ 0.15 EZSI with a negligible phase difference. We have mentioned 
earlier that also the measurement of C^n is hinting at a significant strong phase; leaving this phase as a 
free parameter the allowed range of a from the analysis of the time-dependent CP asymmetry is in good 
agreement with the SM-based indirect estimates of the same. 

The presence of significant strong phase differences in B — > tttt and/or B — > Ktt decays implies that 
the strong interaction effects in these and related decays are not quantitatively described by perturbative 
methods such as QCD factorization. It is, however, less clear if a value of tew ~ 0.3 — 0.5 can also 
be attained in an improved theoretical framework within the SM. So, the current B — ► Ktt data are 
somewhat puzzling. However, judging from the difference R c — R n = 0.41 ± 0.16, and the estimates 
giving R c — R n ~ 0.1, this appears to be about a 2a problem. In view of the fact that current data on 
R c and R n are not quite understood, it is advisable not to use these ratios to constrain 7, as both R c 
and R n involve the poorly understood electroweak penguin contribution. 

To constrain 7 from B — > Ktt data , the so-called mixed ratio Rq, defined below, and advocated some 
time ago by Fleischer and Mannel^SSis potentially useful. The amplitude for the process B + — > K°tt + 
is written in 1|184|) . Using isospin symmetry, the decay amplitude for B° — > K + tt~ can be written as 

A(B° -> K+tt-) = \P'\e lS - |T'|e 17 . (195) 

There is also a small color-suppressed electroweak penguin contribution which can be safely neglected. 
Denoting r = \T'\/\P'\, one has the following relations: 

f (K ± TT T ) 

R ° = f (iW) = 1 ~ 2rCOSSC ° 8 ^ + ^ ' 
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AcHK + n~) ee r^+j+r^-j = -2rsin5sin 7 /Jio, (196) 



which are both functions of r, S and 7. Fleischer and Mannel l zaz l have shown that Rq > sin 2 7 for any 
r and S. So, if Rq < 1, one has an interesting bound on 7. The current value Rq = 0.898 ± 0.071 is 
consistent with one at about 90% C.L., and hence no useful bound emerges on 7. However, if Acp(K + tt~) 
is measured precisely, then the two equations for Rq and RqAcp(K + t: ~) can be used to eliminate S and, in 
principle, a useful constraint on 7 emerges. Current measurements yield Acp{K + n~) = — 0.095 ± 0.029, 
which at 3cr is probably the only significant direct CP-violation observed so far in £?-decays. 

To extract a value of 7 from these measurements (or to put a useful bound), one has to assume a 
value for r, or extract it from data u nder some assumptions. Using arguments based on factorization and 
SU(3)-breaking, Gronau and RosnerESS estimate r from the average of the CLEO, BELLE and BABAR 
data, getting r = 0.142+q'q 2 ^ and a bound 7 < 80° at la. However, there is no bound on 7 at 95% C.L. 
to be compared with the corresponding indirect bounds from unitarity. More precise measurements of 
Ro and Acp(K + tt~) are required to get useful constraints on 7. 

In conclusion, CP asymmetry has been observed in B — ► DK decays using the Dalitz distributions 
at 2.4cr level. Likewise, direct CP asymmetry is seen at 3 a level in Acp{K + it~). However, the current 
significance of the CP asymmetry and the model-dependence of the resonant structure in the decays of 
the .D-meson in the former, and imprecise knowledge of Rq (as well as of Acp(K + ir~)) in the latter, 
hinder at present in drawing quantitative and model-independent conclusions on 7. This situation may 
change with more precise measurements of the various ratios and CP asymmetries in the B — > Kir and 
B — > DK decays. Both the KEK and SLAC £>-factories are now collecting data at record luminosities 
and we trust that improved determinations of 7 (or (^3) will not take very long to come. 



8 Summary and Concluding Remarks 

We have reviewed the salient features of the CKM phenomenology and _B-meson physics, with emphasis 
on new experimental results and related theoretical developments. The data discussed here are spread 
over an energy scale of over three orders of magnitude, ranging from muon decay, determining Gf , to 
the top-quark decays, determining \Vtb\, and have been obtained from diverse experimental facilities. 
Their interpretation has required the development of a number of theoretical tools, with the Lattice- 
QCD, QCD sum rules, chiral perturbation theory, and heavy quark effective theory at the forefront. We 
reviewed representative applications of each of them. Progress in computational technology has enabled 
a quantitative determination of all the CKM matrix elements. While the precision on some of them can 
be greatly improved, all currently available measurements are compatible with the assumption that the 
CKM matrix is the only source of flavour changing transitions in the hadronic sector. In fact, there is 
currently no compelling experimental evidence suggesting deviations from the CKM theory. 

However, there are some aspects of the data which are puzzling and deserve further research. Those 
under current experimental scrutiny are summarized below. 

• Test of unitarity in the first row of the CKM matrix yields Ai = (3.3 ± 1.3) x 10 -3 , which is 2.5 
standard deviations away from zero. Further experimental and theoretical work, yielding robust 
evaluations of the low energy constants of chiral perturbation theory, will have an impact on this 
issue. Precise measurements of decays are being done at DA<E>NE and elsewhere. They will 
yield a determination of \V US \ at significantly better than a per cent level. Analysis of r-decays from 
the B-meson factories will also help. Resolution of the current inconsistencies in the determination 
of qaI 9v m polarized neutron /3-decay experiments will improve the precision on \V u d\ from this 
method. Together, they will enable a precise determination of A^. 

• Experiments at LEP have measured the decays W — > q'q(g), enabling a quantitative test of the 
unitarity involving the first two rows of the CKM matrix. The result Wij\ 2 ~ 2 = 0.039 ± 0.025 is 
consistent with being zero at 1.6 standard deviations. Experiments at CLEO-C and BES-III, but also 
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the B- factory experiments BABAR and BELLE, will measure the matrix elements |V^ S | and \ V c d\ at 
about 1% accuracy, allowing an improved test of the CKM unitarity in the second row. Progress in 
Lattice-QCD technology will be required to have precise knowledge of the D — ► (K, K* ,tt, p) form 
factors. 

• The difference S,pK s — Sj/^K s involving the time-dependent CP asymmetries in the decays B — > <j)Ks 
and B — > J/ipKg, which vanishes in the first approximation in the SM, is currently found to deviate 
from zero. With S^Ks — —0.14 ± 0.33 (not including the scale factor) and S^,k s = —0.14 ± 0.69 
(including the scale factor), the BELLE and BABAR measurements are 1.3 (2.7) standard deviations 
away from Sji^Ks ~ 0.736±0.049 with (without) the scale factor. Including all the b — > sss penguin- 
dominated final states measured so far gives sin 2(3 — (sin2/3 c ff) = 0.50 ± 0.25, with the scale factor, 
which differs from by 2 standard deviations. The presence of a large scale factor (typically 2) in 
the determination of S$k s implies that improved measurements of this (and related) quantities are 
needed to settle the present inconsistency. They will be undertaken at the current and planned B 
factories, and also at LHC-B and B-TeV. 

• Measurements of the ratios R c and R n involving the B — > K~k decays hint at the electroweak penguin 
contributions in these decays at significantly larger strength than their estimates in the SM. However, 
judged from the current measurements, R c — R n = 0.41 ± 0.16, and theoretical estimates yielding 
R c — R n ~ 0.1 based on the default value tew =— 0.15, the mismatch with the SM has a significance 
of about 2 standard deviations. Improved measurements at the B factories and theoretical progress 
in understanding non-leptonic i?-meson decays will clarify the current puzzle. 

• Enhanced electroweak penguin contributions would also influence scmileptonic rare B- and A-decays, 
providing important consistency checks. First round of experiments of the electroweak penguins in 
B — > (K, K* , X s )£ + £ ~ have been reported by the BABAR and BELLE collaborations. Current data 
are summarized ES2i n Table[7|together with the SM-based estimates I ^ ^ of t he same. A comparison 
shows that the experimental measurements are well accounted for in the SM. As the exclusive 
decays have l arger theoretical errors due to the uncertain form factors for which the QCD sum 
rule estimates I I have been used, we quantify a possible mismatch using the theoretically cleaner 
in clusi ve decay B — > X s £ + £~ . Noting that the average of the BELLE and BABAR measurements 



is nMs(B -> X s £+£-) = (6.2 ± l-l±i;l) x 10- 6 , the difference B(X s £+£-) cxp - B(X s £+£-) S m = 



(2.0±2.0) x 10~ 6 amounts to 1 standard deviation. Note that the errors are dominantly experimental. 
From this we infer that there is no evidence of any abnormal electroweak penguin contribution in 
the reliably calculable semileptonic rare -B-decays. Experiments at the B factories and the hadron 
colliders will greatly improve the precision on the decays I! (K, K* , X s )£ + £~ , measuring also 



various distributions sensitive to physics beyond the SM 




Table 7. B — > K^£ + £ and B — > X s £ + £ branching ratios in current experiments ESU and comparison with the SM- 
estimates EHU 



Mode 



BELLE 



BABAR 



Theory (SM) 



B(B 
B(B 
B(B 



K£+£-) (x 10~ 7 ) 
K*£+£-) (x 10- 7 ) 
X s £+£-) (x 10- 6 ) 




6.3 ± 1.6 




± 0.6 



± 1.1 



+ 1.8 
'-1.5 



11.9 ± 3.9 



4.2 ± 0.7 



3.5 ± 1.2 
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From the theoretical point of view, it is very likely that the current deviations from the CKM unitarity 
involving the first row will not stand the force of improved measurements. The anomalies in the penguin- 
dominated £?-decays are also likely to find an experimental resolution, though mapping out the QCD and 
electroweak penguins in i?-decay experiments is crucial in reaching a definitive conclusion. 

Finally, it should be underlined that a number of benchmark measurements in the B- and if-meson 
sectors still remain to be done. On the list of the future experimental milestones are the following: 

• Precise determinations of the weak phases a (or (j> 2 ) and 7 (or 3 ). 

• Measurement of the - B® mass difference AMb s 

• Measurement of the branching ratio B(B® — ► p + p~) 

• Precise measurements of the dilepton invariant mass spectra in B — > (X s , K, K*)£ + £" and the 
forward-backward asymmetries in the decays B — > (X a , K*)£ + £~ . 

• Measurements of the CKM-suppressed radiative and semileptonic rare decays b — > dry and b — > di + £~ 
in the inclusive modes, and some exclusive decays such as B — > (p, iv)-f and B — > (ir, p, u)£ + £~ . 

• Measurements of the rare decays B — ► (X s , K, K*)vv and K — > nvv. 

• Last, but not least, is the challenging measurement of arg(AM Ss ), also called £7, having a value 
57 = — A 2 ?] ~ —2° in the SM. While it appears to be a formidable task to attain the experimental 
sensitivity to probe the SM in 67, searches for beyond-the-SM physics will be undertaken at the 
hadron collider experiments through the CP asymmetry „4cp(-B;s — » J/i>^>). 

There measurements will test the CKM theory in not so well charted sectors where new physics may find 
it easier to intervene. 

It is time to reminisce. Retrospectively, some forty years ago, Cabibbo rotation solved the problem 
of the apparent non- universality of the Fermi weak interactions. The GIM mechanism and the KM 
proposal were crucial steps in the quest of understanding the FCNC processes and CP violation in the 
framework of universal weak interactions. These theoretical developments brought in their wake an entire 
new world of flavour physics. In the meanwhile, all the building blocks predicted by these theories are in 
place. Thanks to dedicated experiments and sustained progress in theory, the field of flavour physics has 
developed into a precision science. All current measurements within errors are compatible with the CKM 
theory. While there is every reason to celebrate and rejoice this great synthesis, this does not necessarily 
imply that we have reached the end of the great saga of discoveries having their roots in flavour physics. 
Looking to the future, one conclusion can already be drawn: As experiments and theoretical techniques 
become more refined, new aspects of physics will come under experimental scrutiny and the known ones 
will be measured with unprecedented precision. It is conceivable that also at the end of the next round 
of experiments, CKM theory will continue to prevail. It is, however, also conceivable that a consistent 
description of experiments in flavour physics may require the intervention of new particles and forces. 
Perhaps, the current experimental deviations from the CKM theory, while not statistically significant, 
are shadows being cast by the coming events. If we are lucky, the hints in the data with renewed efforts 
could turn into irrefutable solid evidence of new physics, or perhaps some of the crucial experiments listed 
above would force us to seek for explanations which go beyond the CKM theory. This remains to be seen. 
Only future experiments can tell if we have reached the shores of a new world or whether the new shores 
that we have reached are still governed by the standard model physics. In any case, the outcome of the 
ongoing and planned experiments in flavour physics will be a vastly improved knowledge about the laws 
of Natute. 
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